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Various  lightning  simulation  test  techniques  were 
conducted  on  the  Lightning  Test  Cylinder  developed  by  the 
Air  Force  Wright  Aeronautical  Laboratories/Atmospheric 
Electricity  Hazards  Group  ( AFWAL/FIESL )  and  the  Air  Force 
Institute  of  Technology  (AFIT/ENG).  The  experimental  tests 
on  the  cylinder,  which  further  investigated  the  assessment 
of  lightning  simulation  techniques  conducted  by  Butters  et 
al.,  included  swept  frequency  continuous  wave  (SFCW), 
current  pulse,  and  shock-excitation.  Designed  to  model  the 
fuselage  of  an  aircraft,  the  aluminum  test  cylinder  is  over 
ten  meters  long  with  a  one  meter  diameter.  To  test  the 
effects  of  various  aircraft  construction  materials,  the 
cylinder  was  constructed  with  an  aperture  where  various 
composite  and  metal  panels  can  be  mounted.  The  research 
involved  determination  of  the  electric  field  and  magnetic 
field  response  transfer  functions  for  each  simulation  test 
technique.  With  these  transfer  functions,  analysis  and 
comparison  of  the  external  and  internal  field  responses 
between  the  SFCW,  current  pulse,  and  shock-excitation  tests 
we^e  made.  A  major  portion  of  the  research  was  to  examine 
the  validity  of  the  linear  model  for  the  current  pulse 
simulation  technique.  In  this  investigation,  transfer 
functions  were  derived  for  various  current  pulse  waveforms. 
The  current  waveforms  injected  into  the  test  cylinder 
included  a  20  kA  unipolar,  double -exponent ia 1  pulse  and  two 
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oscillatory  waveforms  with  peak  amplitudes  of  20  kA  and  100 
kA.  The  research  effort  also  involved  investigation  of  the 
induced  E-dot  (dE/dt)  and  V-dot  (dV/dt)  transients  of  the 
shock-excitation  simulation  test  for  various  composite 
structures.  The  results  of  this  experimental  effort  on  the 
Lightning  Test  Cylinder  were  then  correlated  to  CV-580 
airborne  lightning  measurements.  Differences  and 
similarities  between  the  measured  results  are  presented. 


EXPERIMENTAL  COMPARISON  OF  LIGHTNING  SIMULATION 
TECHNIQUES  TO  CV-580  AIRBORNE  LIGHTNING  STRIKE  MEASURMENTS 


This  thesis  compares  the  Current  Pulse  method  with  the 
Shock-Excitation  method  for  the  simulation  of  the  airborne 
aircraft/lightning  interaction  event.  It  also  provides  a 
comparison  of  these  two  test  techniques  with  Swept  Frequency 
Continuous  Wave  results  and  CV-580  aircraft  in-flight 
measured  lightning  strike  data.  The  major  tasks 
accomplished  during  this  thesis  effort  include: 

a.  Conducting  current  pulse  simulation  tests  on  the 
Lightning  Test  Cylinder  using  unipolar  and  oscillatory 
current  waveforms. 

b.  Conducting  shock-excitation  tests  on  the  Lightning 
Test  Cylinder. 

c.  Developing  transfer  functions  of  magnetic  and 
electric  field  responses  on  the  Lightning  Test  Cylinder 
for  the  current  pulse,  shock-excitation,  and  swept 
frequency  continuous  wave  (SFCW)  simulation  test 
results . 

d.  Comparing  the  temporal  and  spectral  results  of  the 
current  pulse  and  shock -exc i tat i on  simulation  tests  to 


CV-580  aircraft  lightning  strike  data. 
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e.  Computing  the  threat  level  responses  by  folding  the 
threat  current  waveform  into  the  transfer  functions  that 
were  derived  from  the  lightning  simulation  tests. 

Background 

The  interest  in  the  electromagnetic  interaction  of 
lightning  with  aircraft  is  growing  due  to  the  increasing  use 
of  composite  materials,  sensitive  semiconductor  systems,  and 
fly-by-wire  technology,  particularly  in  military  aircraft. 
These  new  aspects  in  aircraft  design  raise  concern  in  the 
susceptibility  and  vulnerability  of  aircraft  to  lightning 
strikes.  Lightning  can  pose  a  serious  hazard  to  the  safety 
and  operation  of  aircraft  with  unprotected  avionic, 
electronic,  and  fuel  systems.  Therefore,  aircraft  design 
requirements  must  include  the  characterization  and 
validation  testing  of  aircraft  protection  from  the  lightning 
threat.  The  test  methods  employed  to  validate  lightning 
protection  must  reproduce  those  significant  electrical 
parameters  required  to  simulate  lightning  effects  on 
aircraft.  Two  common  test  methods  have  been  developed  to 
produce  these  effects;  the  current  pulse  (direct  injection) 
method  and  the  shock-excitation  method.  The  current  pulse 
method  is  most  often  used  by  the  Air  Force  and  the  aerospace 
industry  for  lightning  protection  qualification  testing. 

Some  believe  that  this  method  does  not  produce  the  fast 


changing  electric  field  changes  required  to  fully  simulate 


the  aircraft/lightning  interaction  event.  Although  the 
differences  in  these  two  test  methods  have  been 


investigated,  they  have  never  been  specifically  evaluated  in 
terms  of  actual  in-flight  measured  lightning  strike  data. 

LiaMalna 

Uman  and  Krider  describe  lightning  as  a  transient, 
high-current,  atmospheric  electrical  discharge  (19:79). 
Lightning,  with  its  flash  of  light  and  associated  clap  of 
thunder,  occurs  when  a  particular  region  of  the  atmosphere 
attains  an  electrical  charge  sufficiently  large  that  the 
associated  electric  fields  cause  an  electrical  breakdown  of 
the  air  (19:79).  The  resulting  electric  discharge  may  occur 
within  a  cloud  (intracloud),  cloud-to-cloud  (intercloud),  or 
cloud-to-ground  (19:79). 

Lightning  and  Aircraft 

Aircraft  in  the  vicinity  of  these  highly  charged 
regions  are  susceptible  to  being  struck  by  lightning  and 
becoming  part  of  the  high  current  channel.  An  article  by 
Rustan  reports  that  during  1984  and  1985,  the  Air  Force 
Wright  Aeronautical  Laboratory,  in  conjunction  with  the  FAA, 
U.S.  Navy,  NASA,  and  Office  National  d'Etudes  et  de 
Recherches  Aerospat iales  (ONERA),  obtained  direct  strike 
data  on  an  instrumented  CV-580  aircraft  in  order  to  quantify 
the  lightning  threat  to  aircraft  at  low  altitudes  (17:2).  A 
brief  discussion  of  the  lightning  characterization  program 


is  presented  in  Appendix  A.  The  maximum  levels  of  the 
measured  parameters  were  (17:2): 


1 ) 

peak  current: 

12  kA 

2  ) 

rate-of-rise 

of  current:  3.8  x 

10 

10  A/s 

3  ) 

magnetic  flux 

density 

der i vat i ve 

:  39  50  Wb/m2  s 

4  ) 

electric  flux 

density 

derivative 

:  33  Coul/m2s 

5  ) 

charge  transfer:  103 

coulombs 

b  ) 

pulse  repetit 

ion  rate: 

10  kHz 

During  these  two  years,  the  CV-580  was  only  involved  in  two 
confirmed  cloud-to-ground  lightning  strikes.  Cloud-to- 
ground  strikes  are  believed  to  produce  more  severe  lightning 
effects  than  intracloud  and  intercloud  lightning  strikes. 
Tower  measurements  have  shown  cloud-to-ground  lightning 
parameters  to  be  more  severe  than  those  measured  during  the 
in-flight  strikes  to  the  CV-580  aircraft  (11).  Hebert 
reviewed  the  results  of  three  in-flight  lightning 
characterization  programs  and  found  that  the  rate  at  which 
aircraft  are  struck  by  lightning  decreases  at  lower 
altitudes,  but  the  intensity  of  the  above  parameters 
increases.  As  altitude  increases,  the  rate  of  intracloud 
and  intercloud  strikes  increases  while  the  intensity  of  the 
strikes  decreases  (8).  This  means  that  in  flying  in  or  near 
an  active  thunderstorm,  an  aircraft  is  more  likely  to  be 
struck  at  higher  altitudes,  but  the  intensity  of  such  a 
strike  is  likely  to  be  less  than  at  lower  altitudes.  The 
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lightning  current  parameters  above,  which  are  present  during 
a  lightning  strike,  result  in  both  direct  and  indirect 
effects  to  the  flying  aircraft. 
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The  direct  effects  are  primarily  produced  by  the  charge 
and  energy  transfer  during  the  lightning  strike,  thereby, 
causing  external  physical  damage  to  the  aircraft  (17:4). 

The  electrical  sparkings  that  can  be  produced  during  a 
lightning  attachment,  particularly  inside  the  fuel  tank  of 
an  aircraft,  are  also  considered  direct  effects  (17:4).  An 
article  by  Plumer  and  Robb  reports  that  typical  direct 
effects  on  metallic  structures  include  melting  and 
burnthrough,  pitting  at  structural  interfaces,  resistive 
heating,  magnetic-force  effects,  shock  wave  and  overpressure 
damage,  and  arcing  across  bonds,  hinges,  and  joints 
(14:159-163).  On  nonmetallic  structures,  the  direct  effects 
of  lightning  strikes  include  puncturing  of  fiberglass  or 
Kevlar  reinforced  plastics  and  the  delamination  and  strength 
degradation  of  graphite  composite  structures  (14:163-167). 

The  indirect  effects  of  lightning  on  aircraft  are  the 
induced  electromagnetic  interactions  with  aircraft  wiring 
and  avionics  systems  (14:1585.  This  can  include  inductive 
and  capacitive  coupling  mechanisms.  The  indirect  effects 
are  determined  from  the  relationship  between  the  surface 
currents  and  charge  densities  throughout  the  aircraft 
surface  and  the  induced  transients  on  the  wiring  and 
avionics  systems  (17:4). 
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Perala  and  DuBro  point  out  that  aircraft/lightning 
interaction  events  are  the  result  of  the  high  voltage  phase 
and  high  current  phase  of  the  lightning  strike  (12:3).  The 
high  voltage  phase  occurs  when  the  lightning  channel 
initially  attaches  to  the  aircraft  (12:3).  During  this 
phase,  there  is  a  rapid  rise  and  time  rate-of -change  in  the 
electric  fields  and  a  relatively  small  current  induced  on 
the  aircraft  surface  (12:3).  This  is  followed  by  the  high 
current  phase,  which  occurs  when  the  return  stroke 
propagates  through  the  lightning  channel  (12:3). 

Lightning  and  Simulation 

Because  of  the  complexity  of  the  aircraft/lightning 
interaction  process,  it  is  necessary  to  use  ground  lightning 
simulation  and  testing  to  assess  the  effects  of  lightning  on 
aircraft.  Laboratory  lightning  simulation  does  not  attempt 
to  actually  reproduce  the  entire  natural  lightning  event, 
since  natural  lightning  is  a  highly  varying,  discontinuous, 
and  complex  phenomenon.  Rather,  the  purpose  of  lightning 
simulation  is  to  generate  those  electrical  parameters  which 
reproduce  the  direct  and  indirect  effects  of  lightning  on 
aircraft.  A  report  by  Butters  et  al.  (2)  describes  two 
major  lightning  simulation  techniques  which  have  been 
developed  to  produce  lightning's  most  significant  electrical 
parameters:  current  pulse  and  shock-excitation. 

The  current  pulse  technique  is  performed  by  injecting  a 
low-level  (or  high-level)  current  pulse  waveform  through  the 


-6- 


test  aircraft  (2:1).  In  the  low-level  case,  the  resulting 
induced  transients  are  measured  and  then  linearly  scaled  to 
the  required  threat  levels  (2:1).  When  moderate  threat 
waveforms  (30,000  Amperes)  are  injected,  only  minor  scaling 
is  required.  When  full  threat  pulses  (200,000  Amperes)  are 
injected,  no  scaling  is  required.  It  should  be  noted  that 
there  are  only  two  lightning  simulators  in  the  United  States 
capable  of  injecting  the  full  threat  pulse  into  a  fighter 
sized  aircraft  or  larger:  one  at  Sandia  National 
Laboratories  and  one  developed  by  the  Air  Force  in  its 
Atmospheric  Electricity  Hazards  Protection  Advanced 
Development  Program  (7).  The  majority  of  current  pulse 
lightning  simulation  tests  are  performed  at  reduced  current 
levels  and  the  resulting  transients  are  linearly  scaled  to 
the  full  threat  level.  With  this  technique,  a  current 
impulse  generator  and  a  waveforming  network  are  attached  to 
the  aircraft  test  body  (2:3).  The  current  return  path  is 
then  wired  from  the  aircraft  back  to  the  generator  (2:3). 

The  shock-excitation  technique  is  a  recent  development 
in  lightning  simulation  technology.  This  technique  uses 
spark  gaps  to  isolate  the  aircraft  test  body  from  the 
generator  and  ground  (2:1).  The  shock-excitation  technique 
discharges  a  high  voltage  generator  (a  Marx  bank)  to  break 
down  the  air  between  the  gaps  and  to  provide  both  current 
and  voltage  excitation  (2:1).  The  result  is  reported  to 
reproduce  the  effects  of  several  lightning  strike  events: 


the  rapid  change  of  the  electric  fields,  the  lightning 
attachment,  and  the  return  stroke  (2:1). 

In  comparing  the  two  lightning  simulation  techniques. 
Butters  points  out  that  the  current  pulse  technique  assumes 
dominant  inductive  coupling  mechanisms  and  the  absence  of 
nonlinear  lightning  effects  (2:1).  However,  higher  current 
levels  result  in  nonlinear  effects  and  deviate  from  the 
linear  scaling  analysis  (2:3).  But,  it  is  argued  that  these 
deviations  will  result  in  a  more  conservative  linear 
analysis  of  the  lightning  effects  on  the  aircraft  (2:3). 
Perala  and  DuBro  state  that  the  current  pulse  technique 
simulates  only  the  high  current  phase  of  the  lightning 
strike;  whereas,  the  shock-excitation  technique  simulates 
both  the  high  current  and  high  voltage  phases,  and  also 
addresses  the  analysis  of  the  nonlinear  effects  (12:3,5). 

Riablcro 

The  overall  goal  of  this  thesis  is  to  assess  the 
capabilities  of  the  current  pulse  and  the  shock-excitation 
simulation  methods  to  produce  the  significant  lightning 
electrical  parameters,  which  contribute  to  the  airborne 
aircraft/lightning  interaction  event.  Specifically,  the 
objective  is  to  determine  the  answers  to  the  following 
questions : 

a.  How  do  the  transients  for  the  current  pulse  and 


shock-excitation  simulation  techniques  compare  to  swept 


frequency  continuous  wave  measurements  when  linearly 
scaled  to  the  full  threat  level? 

b.  Is  the  current  pulse  method  alone  sufficient  to 
assess  the  susceptibility/vulnerability  of  an  aircraft 
to  the  lightning  threat? 

c.  Does  the  current  pulse  method  faithfully  reproduce 
all  of  the  significant  effects  recorded  during  the 
in-flight  lightning  strikes  to  the  CV-580? 

d.  What  significant  effects,  if  any,  are  not 
reproduced  by  the  current  pulse  method  and  can  only 
be  produced  by  using  the  shock-excitation  method? 

Scape 

Although  there  are  several  variations  associated  with 
lightning  simulation  techniques,  this  thesis  will  only 
investigate  the  current  pulse  simulation  technique  using 
unipolar  and  oscillatory  waveforms,  the  shock -exc i tat i on 
simulation  technique,  and  the  SFCW  technique.  The 
simulation  tests  are  conducted  on  the  Lightning  Test 
Cylinder.  Also,  there  have  been  several  in-flight  lightning 
character ization  programs  prior  to  the  CV-580  program  during 
1984  and  1985.  The  current  pulse  and  shock-excitation 
simulation  techniques  will  only  be  compared  with  the  data 
measured  and  recorded  by  the  CV-580  lightning 
character i zat ion  program.  The  parameters  that  will  be  used 
for  comparison  in  this  thesis  are  limited  to  the  maximum 
measured  and  computed  levels  of  the  E-field,  H-field,  time 


rate-of -change  of  the  E-field  (dE/dt),  time  rate-of -change 
of  the  H-field  (dH/dt),  and  current.  Frequency  domain 
comparisons  will  only  be  directed  to  spectral  distributions 
of  the  transfer  functions  for  the  airborne  measurements  and 
each  of  the  simulation  techniques. 

Assumptions 

One  key  assumption  necessary  for  this  thesis  is  that 
the  induced  transients  of  the  aircraft/lightning  interaction 
event  can  be  lii^  rly  extrapolated  to  the  full  threat  level. 
Therefore,  the  interaction  event  can  be  modeled  as  a  linear, 
time  invariant  system  with  the  use  of  transfer  functions. 
Another  required  assumption  is  that  a  digital  sample  rate  of 
5  nanoseconds  for  a  10  microsecond  window  is  sufficient  for 
the  data  acquisition  equipment  to  record  the  significant 
interaction  effects  during  the  in-flight  tests  and  the 
lightning  simulation  tests.  Also,  an  assumption  is  made 
that  the  CV-580  data  is  representative  of  a  lightning 
interaction  event  with  an  airborne  aircraft. 

Appg.Qac.ti 

The  approch  of  this  thesis  is  to  experimentally  and 
analytically  determine  which  test  method  most  faithfully 
reproduces  lightning  interaction's  most  significant  effects 
by  making  both  time  domain  and  frequency  domain  comparisons 
of  the  current  pulse  and  shock -exc i tat i on  tests  to  the 


CV-580  airborne  measurements. 


The  first  part  of  the  thesis 


involves  conducting  experimental  tests  on  the  Lightning  Test 
Cylinder  developed  by  the  Air  Force  Wright  Aeronautical 
Laboratories  --  Atmospheric  Electricity  Hazards  Group 
( AFWAL/FIESL )  and  the  Air  Force  Institute  of  Technology 
(AFIT).  The  experimental  tests  on  the  test  cylinder  include 
current  pulse  and  shock-excitation  lightning  simulations  and 
swept  frequency  continuous  wave  measurements.  From  these 
experimental  measurements,  the  response  transfer  functions 
of  the  test  cylinder  from  each  test  are  then  developed. 

With  the  transfer  functions,  analysis  and  comparisons  of  the 
responses  beween  the  SFCW,  current  pulse,  and  shock- 
excitation  tests  to  the  threat  level  lightning  waveform  can 
be  made . 

The  second  part  of  this  thesis  is  comparison  to  the 
CV-580  airborne  lightning  strike  data.  The  transfer 
function  of  the  aircraft's  response  to  an  actual  lightning 
strike  is  compared  to  those  of  the  current  pulse  and 
shock -exc i tat  ion  lightning  simulation  methods.  Lightning 
interaction  effects  are  analyzed  by  comparisons  and 
evaluations  of  transfer  functions  between  induced  transients 
external  and  internal  to  the  test  cylinder. 

Sequence  of  Presentation 

This  thesis  is  developed  in  several  distinct  sections. 
In  Chapter  2,  the  elementary  theory  necessary  to  understand 


the  basic  interaction  process  of  lightning  with  aircraft  and 


r 


the  lightning  simulation  techniques  is  presented.  The 
experimental  procedures  and  data  acquisition  of  the 
lightning  simulation  measurements  are  discussed  in  Chapter 
3.  Chapter  4  presents  comparisons  and  analysis  of  the 
experimental  results  with  the  airborne  lightning  strike 
data.  Conclusions  and  recommendations  are  discussed  in 
Chapter  5. 
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An  understanding  of  the  basic  interactions  of  lightning 


with  aircraft  and  the  simulation  techniques  to  reproduce 
lightning's  significant  effects  is  necessary  for  the 
development  of  this  thesis.  This  chapter  will  present  and 
review  the  lightning/aircraft  electromagnetic  interaction 
process,  lightning  simulation  requirements,  the  current 
pulse  technique,  the  shock-excitation  technique,  and 
frequency  domain  analysis  methods. 

Lightning  Interaction  with  Aircraft 

The  interest  in  the  electromagnetic  interaction  of 
lightning  with  aircraft  is  increasing  due  to  the  growing  use 
of  composite  materials,  sensitive  semiconductor  systems,  and 
fly-by-wire  technology.  This  section  will  review  the 
fundamentals  of  the  interaction  process.  An  extensive 
review  of  state  of  the  art  lightning  interaction  modeling 
has  been  written  by  Perala,  Rudolph,  and  Eriksen  (13). 

The  1 ightning/a ircra f t  e lectromagnet ic  interaction 
process  is  often  considered  as  three  separate  processes: 
external  interaction,  internal  interaction,  and  internal 
propagation.  Although  these  three  processes  are  analyzed 
individually,  theoretically  they  are  not  separate  (13:173). 
The  interaction  process  can  be  analyzed  separately,  however, 
when  there  is  very  little  mutual  coupling  between  the 
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Figure  1.  Lightning/Aircraft  Interaction  Process  (13) 


The  first  part  of  the  process,  the  external 
interaction,  is  initiated  when  the  lightning  channel  induces 
a  surface  current  density  and  a  charge  density  which  results 
from  the  flow  of  electrical  charge  on  the  aircraft  surface. 
The  surface  current  and  charge  densities  respectively 
correspond  to  the  tangential  magnetic  fields  and  normal 
electric  fields  on  the  surface  of  the  aircraft  (13:173). 


Internal  interaction  is  the  second  part  of  the 
lightning/aircraft  interaction  process.  It  is  the 


penetration  of  lightning  energy  into  the  interior  of  the 
aircraft  due  to  apertures,  antennas,  exposed  cables  and 
conductors,  panel  seams  and  joints,  and  skin  diffusion  of 
the  aircraft  (13:173). 

The  final  part  of  the  interaction  process,  internal 
propagation,  follows  once  the  lightning  energy  penetrates 
through  the  shell  of  the  aircraft.  In  general,  the 
aircraft's  internal  wiring,  metallic  fluidlines,  and  control 
cables  act  as  transmission  lines  through  which  the  lightning 
energy  effectively  propagates,  ultimately  to  critical 
electronic  systems  (13:173). 

Lightning  Simulation  and , Testing. 

Ground  lightning  simulation  and  testing  are  used  to 
assess  the  effects  of  lightning  on  aircraft  and  to  validate 
lightning  protection  for  aircraft.  The  goal  of  lightning 
simulation  is  not  to  reproduce  the  entire  natural  lightning 
event,  but  to  accurately  simulate  or  reproduce  the 
significant  direct  and  indirect  effects  of  lightning  on 
aircraft.  The  simulation  of  these  direct  and  indirect 
effects  are  often  based  on  the  production  of  four  of 
lightning's  electrical  parameters  (6:2): 

1)  the  current  peak  amplitude  (Imax) 

2)  the  maximum  current  rate  of  change  (dl/dt) 

3)  the  charge  transfer  (  /  idt) 


4)  the  action  integral  (  /i2dt) 

The  action  integral  (  / i2dt)  represents  the  lightning 
pulse's  ability  to  dissipate  or  deliver  energy  to  the 
aircraft  or  test  body  (3:209).  In  general,  lightning 
simulation  and  testing  programs  incorporate  statistical 
combinations  of  the  four  above  parameters  in  order  to 
characterize  the  lightning  interaction  event  (6:2). 

Criteria  for  military  lightning  simulation  testing  is 
documented  in  MIL-STD-1757A,  which  presents  the  waveforms  to 
be  used  in  laboratory  lightning  simulation  tests  and  for 
lightning  protection  qualification  of  military  aircraft. 
MIL-STD-1757A  specifies  the  following  test  waveforms 
(3:210) : 

1)  Voltage  Waveform  A:  for  punch-through/f lashover 
testing;  a  1000-kV/ps  ramp  until  breakdown. 

2)  Voltage  Waveform  B:  dielectric  testing  and  fuel 
ignition  testing;  1.2  jus  time-to-peak  and  50  ps 
t ime-to-decay  to  50  percent. 

3)  Voltage  Waveform  C:  attach-point  testing;  2  ps 
time-to-peak . 

4)  Voltage  Waveform  D:  attach-point  testing;  50  ps 
time-to-peak . 

5)  Current  Waveform  A:  initial  stroke  component. 

6)  Current  Waveform  B:  intermediate  current  component 

7)  Current  Waveform  C:  continuing  current  component. 

8)  Current  Waveform  D:  restrike  component. 

9)  Current  Waveform  E*  magnetic-coupling  testing. 


Typically,  these  waveforms  are  utilized  and  applied  to  test 
aircraft  by  appropr iately  time  sequencing  the  various 
waveform  components  (3:210).  The  test  waveforms  are  not 
designed  to  replicate  natural  lightning  itself,  but  to 
simulate  the  direct  and  indirect  effects  caused  by  the 
lightning  event  (3:211).  Note  that  while  many  lightning 
simulation  tests  which  employ  the  current  pulse  (direct 
injection)  method  concentrate  primarily  on  the  current 
waveform  components  presented  above,  the  military  standard 
includes  voltage  waveforms  as  well. 

Current  Pulse  Technique 

The  current  pulse  technique,  also  known  as  the 
lightning  transient  analysis  ( LTA)  technique,  was  originally 
developed  by  the  General  Electric  High  Voltage  Laboratory 
(15:15).  It  is  performed  by  injecting  a  current  pulse  — 
low-level,  moderate-level,  or  full  threat-level  --  through 
the  structure  of  the  test  aircraft  (2:1).  With  this 
technique,  a  current  impulse  generator  (a  capacitor  bank) 
and  a  resistive  waveshaping  network  are  attached  directly  to 
the  aircraft  test  body  at  an  entry  point  (2:3;  21:450).  The 
current  return  path  for  the  setup  is  then  hard-wired  from 
the  exit  point  on  the  aircraft  back  to  the  generator  (2:3). 
The  test  configuration  is  often  modelled  as  a  simple  series 
RLC  circuit  (21:450).  Analytically,  the  inductance  of  the 
configuration  is  assumed  to  be  determined  by  a  combination 


of  the  current  return  path  and  the  test  aircraft  (21:450). 
The  waveshaping  network  and  the  current  generator  determine 
the  resistance  and  the  capacitance  of  the  system 
respectively  (21:450).  An  example  of  the  current  pulse 
setup  --  a  test  on  the  F-16  by  AFWAL  --  and  its  equivalent 
circuit  are  illustrated,  respectively,  in  Figures  2  and  3. 
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Figure  2.  Current  Pulse  Test  Configuration  (21:456 


Figure  3.  Current  Pulse  Equivalent  Circuit  (15:16) 


The  parameters  of  the  configuration  are  chosen  such 
that  a  unipolar,  double  exponential,  impulse  current 
waveform  is  injected  through  the  aircraft  structure 
(21:450).  The  equations  describing  the  injected  current 
waveform  are  given  by  (21:450) 


Ik  -  alt  -Ik  -  alt 

e  -  e 


_  R 


k  =  a 


where 


i ( t )  =  time  varying  current  (amperes) 
V  =  capacitor  voltage  (volts) 

L  =  circuit  impedance  (henries) 

R  =  circuit  resistance  (ohms) 


C  =  circuit  capacitance  (farads) 


The  resulting  normalized  waveform,  marked  by  the  rise  time 
(tr)  and  decay  time  (tt),  is  illustrated  in  Figure  4. 

Most  often  the  current  pulse  technique  is  used  to 
inject  a  low  threat-level  or  moderate  threat-level  current 
waveform  from  which  the  resulting  induced  transients  are 
linearly  extrapolated  to  full  threat  levels,  typically  a 
current  peak  of  200  kA  (21:449).  In  order  to  simulate  the 


Figure  4.  Unipolar  Test  Waveform  (21:456). 


lightning/aircraft  interaction  event,  MIL-STD-1757A  and  data 
from  statistical  lightning  character i zat i ons  reguire  the 
following  criteria  for  the  current  pulse  waveform  (21:451): 

Ip  =  30  kA  (moderate  threat-level  peak  current) 
tr  =  2  /js  +.  20%  (10%  to  90%  rise  time) 

t(  =  50  )is  +  20%  (50%  tail  time) 

The  resistive,  capacitive,  and  inductive  components  of  the 
test  configuration  are  selected  to  maintain  the  peak  current 
IP  while  keeping  the  waveform  within  the  required  tolerances 
(21:452).  Computer  numerical  techniques  are  used  to 
optimize  the  test  configuration  by  determining  the  required 


values  of  the  components  (21:452).  Using  the  Kirchoff 
voltage  equation,  the  test  configuration  equivalent  circuit 
produces  the  following  series  RLC  equation  at  time  t  =  0+ 
(21:452) : 

V  =  (l/C,/i(t)dt  +  R  i ( t )  +  L  di(t)/dt  (4) 

whose  time  derivative  is  given  as: 

0  =  (1/C)  i ( t )  +  R  di(t)/dt  +  L  d2 i ( t ) /dt 2  (5) 

For  full  scale  lightning  tests  on  the  F-16  aircraft, 
the  current  pulse  test  setup  at  AFWAL ' s  Flight  Dynamics 
Laboratory  uses  two  banks  of  capacitors  with  a  total 
capacitance  of  8.4  microfarads  (21:453).  The  configuration 
is  improved  by  utilizing  a  pneumatic  spark  gap  system  that 
closes  until  sparkover  between  the  two  banks  to  produce  a 
cleaner  current  waveform  (15:19;  21:454).  The  current  is 
then  damped  by  a  resistive  network  of  3  ohms  at  the 
generator  output  (21:453).  Another  improvement  of  the 
current  pulse  technique  at  the  Air  Force  Flight  Dynamics 
Laboratory  is  the  use  of  multiple  return  paths  (15:19; 
21:450-451).  The  multiple  paths  more  evenly  distribute  the 
field  around  the  test  aircraft  and  reduce  the  inductance  of 
the  system  to  obtain  a  higher  peak  current  (21:450-451).  By 
utilizing  four  aluminum  return  paths,  the  overall  inductance 


of  the  configuration  with  the  F-16  aircraft  is  approximately 
3.5  microhenries  (21:454).  In  measuring  the  induced 
transients  from  lightning/aircraf t  coupling,  the  current 
pulse  simulation  setup  at  the  Flight  Dynamics  Laboratory 
incorporates  a  transient  digitizer  to  provide  spectral  data 
on  the  phase  and  amplitude  of  the  driving  current  and 
response  (15:19).  From  this  information,  the  actual 
transfer  function  of  the  test  configuration  is  then 
determined  and  appropriately  identified  in  the  analysis  of 
the  results  (15:19). 

In  the  current  pulse  lightning  simulation  technique, 
the  resulting  measured  induced  transients  from  the  injected 
low-level,  or  moderate-level,  current  pulse  are  linearly 
scaled  to  full  threat  levels  where  res istively-coupled 
transients  are  extrapolated  from  the  peak  current  and 
inductively-coupled  transients  are  extrapolated  from  the 
current  time  rate-of -change .  The  current  pulse  technique 
assumes  that  linear  extrapolation  of  the  induced  transients 
is  valid  if  the  injected  current  waveform  is  identically 
scaled  to  the  full  threat  lightning  waveform  (2:1,3). 

From  analytical  and  experimental  investigations, 
McCormick,  Maxwell,  and  Finch  report  on  the  validity  of  the 
current  pulse  (LTA)  technique  (10).  In  the  research  program 
sponsored  by  the  Air  Force  Flight  Dynamics  Laboratory, 
lightning  simulation  test  (LST)  data  is  used  to  derive  the 


transfer  function  of  the  test  configuration  (10).  Using  the 


time  and  frequency  domain  inputs  derived  from  the  LST  data, 
it  was  concluded  that  the  LTA  technique  is  a  valid  method 
for  assessing  the  susceptibility  and  vulnerability  of 
aircraft  to  the  effects  and  hazards  of  lightning  (10:142). 
The  results  also  determined  that  the  linearity  of  the 
current  pulse  configuration  is  a  valid  and  justified 
assumption,  and  that  linear  extrapolation  to  the  full  threat 
lightning  level  is  an  appropriate  analysis  (10:134,142). 

Shock -Excitation  Technique 

The  shock-excitation  technique,  developed  by  the 
McDonnell  Aircraft  Company,  differs  from  the  current  pulse 
technique  by  using  spark  gaps  to  electrically  isolate  the 
aircraft  test  body  from  the  generator  and  ground  (2:1; 
15:23).  The  shock-excitation  setup  incorporates  these  large 
input  and  output  spark  gaps  on  the  test  aircraft  to  produce 
both  charging  transients  and  discharging  transients.  A  high 
voltage  Marx  generator  is  used  to  break  down  the  air  between 
the  spark  gaps  and  to  provide  both  current  and  voltage 
excitation  to  the  aircraft  (2:1).  A  Marx  bank  is  a  simple 
pulsed  power  generator  in  which  N  capacitors  are  charged  in 
parallel  to  a  given  voltage,  V0 .  The  capacitor  bank  is  then 
discharged  in  series  to  achieve  very  high  voltages;  ideally, 
N  x  VQ  volts.  The  arc  breakdown  of  air  occurs  at 
approximately  10  kilovolts  per  centimeter.  Therefore,  an 
airgap  of  1  meter  requires  approximately  one  million  volts 


to  break  down. 


In  the  sequence  of  events  during  the  shock -exc i tat i on 


simulation,  the  discharging  Marx  generator  produces  a 
rapidly  changing  E-field  between  the  input  electrode  and  the 
test  body.  After  a  few  microseconds,  an  arc  will  be 
established  which  will  produce  a  rapid  charging  of  the  body. 
A  second  arc  will  then  be  established  from  the  body  to 
ground  --  rapidly  discharging  the  vehicle.  The  two  arcs 
will  complete  the  circuit  of  the  test  setup  and  allow  the 
Marx  generator  to  discharge  a  high  current  through  the 
vehicle.  Because  of  the  finite  amount  of  time  required  for 
the  long  arcs  to  be  established,  various  transient 
mechanisms  are  simulated.  Reportedly,  the  different 
mechanisms  simulated  with  the  shock-excitation  technique 
are:  (1)  nearby  lightning  from  the  rapidly  changing 

E-field,  (2)  stepped-leader  attachment  from  the  charging  of 
the  aircraft,  and  (3)  return  stroke  from  the  discharging  of 
the  test  vehicle  (2:1).  The  McDonnell  Aircraft  Company 
shock -exc i tat i on  test  setup  is  illustrated  in  Figure  5. 

The  shock-excitation  configuration  and  the  resulting 
coupling  to  internal  wiring  are  modeled  analytically  as  two 
coupled  transmission  line  circuits,  as  shown  in  Figure  6, 
where  one  circuit  represents  the  aircraft  test  body  and  the 
other  represents  the  internal  wiring  of  the  aircraft 
(2:12-13).  The  generator  and  the  waveshaping  network  are 
modeled  as  lumped  parameters  (2:12).  A  time-varying  arc 
resistance  is  used  to  model  the  spark  gaps  in  the  test 
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configuration  (2:12).  The  capacitive  and  inductive  coupling 
mechanisms  are  modeled,  respectively,  as  a  parallel  current 
source  Is  and  a  series  voltage  source  Vs  in  the  equivalent 
circuit  of  the  internal  wiring  (2:13). 

The  shock-excitation  lightning  simulation  can  be 
divided  into  a  charging  phase  and  a  discharging  phase  (2). 
The  charging  phase  in  essence  simulates  the  stepped-leader 
attachment  of  the  lightning  strike,  and  the  discharging 
phase  simulates  the  return  stroke  (2:37).  The  charging 
phase  is  simulated  by  isolating  the  aircraft  test  body  from 
ground  potential  and  then  switching  the  high-voltage 
generator  output  to  the  test  body  (2:37).  The  high  voltage 
potential  of  the  generator  induces  the  flow  of  displacement 
currents  on  the  surface  of  the  test  body,  thus,  quickly 
charging  the  aircraft  structure  to  a  high  potential  (2:37). 
The  resulting  currents  and  voltages  on  the  test  body  during 
the  charging  phase  is  modeled  from  a  series  RLC  circuit 
(2:15,38).  The  current  and  voltage  waveforms  are  given 
respectively  as  (2:15,38): 

i(t)  =  (V/wL)  e'at  sin(wt)  (6) 

and 

V(t)  =  (1/Cc)  J  i  ( t ) dt  (7) 
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where 


a  =  R/2L 

w  =  frequency  (radians/sec) 

V  =  generator  voltage  (volts) 

C c  =  test  body  capacitance  (farads) 

L  =  total  inductance  (henries) 

R  =  series  resistance  (ohms) 

The  time  derivatives  of  the  current  and  voltage  waveforms 
are  given  respectively  as  (1:38): 


di (t)/dt 


(V/wL) 


( w) cos ( wt ) - ( a ) s in ( wt ) 


-  a  t 


e 


(8) 


and 


dV(t)/dt  =  i(t)/Cc 


(9  ) 


Following  the  charging  phase,  the  discharging  phase  of 
the  shock-excitation  technique  simulates  the  return  stroke 
of  the  lighr.ning  strike  (2:51).  During  the  charging  phase, 
the  test  body  is  charged  to  a  high  potential  due  to  high 
electrical  fields  produced  by  the  Marx  generator's  output 
and  due  to  the  output  spark  gap  isolating  the  test  body  from 
ground  potential  (2:51).  The  discharge  phase  occurs  once 
the  output  spark  gap  breaks  down  and  the  test  body  is  then 
quickly  discharged  to  ground  potential  (2:51).  This  rapid 
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discharge  produces  very  high  dE/dt's  of  approximately 


1012  V/m/s  (2:51). 

Frequency  Domain  Analysis 

Frequency  domain,  as  well  as  time  domain,  techniques 
are  used  in  this  thesis  for  the  analysis  of  measurements 
from  the  CV-580  airborne  lightning  program,  the  current 
pulse  tests,  and  the  shock-excitation  tests.  In  the  first 
part  of  the  analysis,  frequency  domain  techniques  are 
applied  to  the  measured  data  to  determine  the  actual 
response  of  both  the  airborne  and  simulated 

lightning/aircraft  interaction  events.  The  measured  sensor 
data  not  only  records  the  electrical  phenomena,  but  also  the 
particular  characteristics  of  the  sensor.  A  more  accurate 
characterization  of  the  electrical  phenomena  is  provided  by 
isolating  and  removing  the  effects  of  the  measurement  system 
from  the  data  (6:5).  First,  frequency  domain  transfer 
functions  of  each  component  in  the  measurement  system  are 
determined  by  a  Hewlett-Packard  3577A  network  analyzer.  The 
components  of  the  measurement  system  include  electric  field, 
magnetic  field,  and  current  sensors  with  fiber  optic  links. 
The  network  analyzer  measures  the  response  of  each 
measurement  system  component  to  a  swept  frequency  continuous 
wave  ( SFCW )  signal,  from  DC  to  the  upper  frequency  limit 
associated  with  the  component  (6:5).  The  ratio  of  the 
response  to  the  CW  signal  at  discrete  frequencies  produces 


the  corresponding  frequency  domain  transfer  function  of  each 
component  in  the  measurement  system  (6:5). 

The  measured  time  domain  transients  due  to  lightning 
interaction  effects  are  also  transformed  into  the  frequency 
domain  by  using  the  Fast  Fourier  Transform.  Dividing  out 
the  transfer  functions  of  the  measurement  system  components 
from  these  measured  transients  produces  the  actual  lightning 
response  in  the  frequency  domain.  An  inverse  of  this 
corrected  frequency  transform  results  in  the  actual  time 
domain  signal  by  removing  the  effects  of  the  measurement 
system  (6:5). 

An  advantage  of  this  procedure  is  that  the  response  of 
a  derivative  field  sensor,  electric  or  magnetic,  is 
integrated  in  the  time  domain,  as  well  as  removing  the 
effects  of  the  measurement  system.  In  the  frequency  domain, 
dividing  out  the  transfer  function  of  the  derivative 
measurement  component  corresponds  to  integration  in  the  time 
domain. 

The  second  part  of  the  frequency  domain  analysis  is  to 
produce  the  transfer  functions  of  the  CV-580  aircraft  and 
the  test  object  at  the  various  sensor  locations.  The 
transfer  functions  that  are  developed  relate  the  response  of 
the  system  to  airborne  or  simulated  lightning  strikes.  This 
is  shown  in  the  following  relationship  (6:5): 


T(«)  =  R(w)/S( 
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where 


T(~)  =  frequency  domain  transfer  function 

R(w)  =  Fourier  transform  of  actual  sensor  response 


S ( ^ )  =  Fourier  transform  of  applied  excitation  source 

The  response  transfer  function  T(oj)  contains  both  amplitude 
and  phase  information.  They  are  produced  for  the  SFCW, 
current  pulse,  shock-excitation,  and  airborne  cases  at  each 
sensor  location.  The  excitation  source,  as  well  as  the 
sensor  response,  may  be  voltage,  current,  electric  field,  or 
magnetic  field. 

Once  the  transfer  function  is  determined  it  is  possible 
to  analytically  determine  what  the  transient  response  at  a 
particular  location  would  be  to  various  threat  waveforms. 
This  is  shown  by  the  following  relationship: 

Rt(«)  =  IT  {<-)  T(w)  (11) 


wher  e 

R  T { — )  =  response  to  threat  waveform 
IT(^)  =  threat  waveform 
T  (  -*-  )  =  transfer  function 

By  taking  the  inverse  Fourier  transform  of  RT(~),  it  is 
possible  to  obtain  the  time  domain  transient  response  to  the 
threat  waveform. 
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It  is  important  to  realize  that  the  frequency  domain 
analysis  using  the  Fourier  transform  assumes  a  linear 
t ime - invar iant  system.  Therefore,  the  key  assumption  in  the 
analysis  is  that  the  lightning  interaction  event  is  a  linear 
process  and  neglects  nonlinearities  such  as  arcing,  corona, 
streamering,  and  surface  tracking.  It  is  argued,  however, 
that  these  mechanisms  result  in  losses  of  electromagnetic 
energy  and  that  an  assumption  of  linearity  will  result  in  a 
more  conservative  analysis  of  the  lightning  interaction 
event  --  predicting  higher  values  than  if  the  nonlinearities 
are  taken  into  cons iderat ion  (2:3).  Actual  tests  at  very 
low  level  SFCW  inputs  of  about  25  watts,  followed  by  current 
pulse  tests  at  20,000  amperes,  showed  that  the  extrapolation 
was  conservative  in  the  majority  of  cases  during  lightning 
simulation  tests  on  an  F-14  fighter  aircraft  (5).  Because 
of  the  exceptions  to  the  case  and  the  importance  of 
validated  lightning  protection,  extreme  care  must  be 
exercised  in  applying  frequency  domain  techniques  for  actual 
lightning  protection  qualification  testing  when  low  level 
currents  are  used.  In  general,  when  applied  to  transient 
reponses  recorded  during  current  pulse  tests,  the  less  the 
extrapolation  necessary,  the  better  (5). 
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3.  Lightning  Simulation  Measurements 

This  chapter  describes  the  experimental  phase  of  this 
thesis.  It  discusses  the  procedures  and  equipment  of  the 
simulation  tests  that  were  conducted  on  the  Lightning  Test 
Cylinder  at  AFWAL/FIESL.  The  lightning  simulation 
techniques  conducted  on  the  Lightning  Test  Cylinder  were 
Swept  Frequency  Continuous  Wave,  Current  Pulse,  and 
Shock-Excitation.  The  resulting  magnetic  fields,  electric 
fields,  and  currents  induced  on  the  cylinder  are  recorded 
for  each  test  technique. 

Lightning  Test  Cylinder 

The  Lightning  Test  Cylinder  is  the  test  bed  that  was 
used  for  the  lightning  simulation  experiments  performed  for 
this  thesis.  The  cylinder  was  fabricated  by  the  Air  Force 
Institute  of  Technology  Fabrication  Shop  out  of  0.060  inch 
thick  sheet  aluminum  using  aircraft  construction  techniques. 
Designed  to  model  the  fuselage  of  an  aircraft,  the  cylinder 
is  over  ten  meters  in  length  with  a  diameter  of  one  meter. 

In  order  to  determine  the  shielding  effectiveness  and  the 
entry  of  electromagnetic  energy  through  various  aircraft 
construction  materials,  the  cylinder  was  built  with  an 
aperture  in  the  center.  The  aperture  is  1.5  meters  long  and 
spans  one-third  of  the  circumference  of  the  cylinder.  A 
variety  of  panels  can  then  be  fastened  over  this  aperture. 


The  panels  used  in  this  thesis  include  one  solid  aluminum 


panel,  which  is  identical  to  the  material  of  the  cylinder, 
and  two  graphite  composite  panels.  One  of  the  composite 
panels  is  composed  of  nickel-coated  carbon  fibers  and  is 
identified  in  this  thesis  as  Composite  #1.  The  other 
composite  panel,  which  is  identified  as  Composite  #2,  is 
also  composed  of  nickel-coated  carbon  fibers  and  has  a 
copper  mesh  imbedded  in  the  surface  of  the  panel.  The  tests 
were  also  conducted  with  an  open  aperture,  using  no  panel. 


A  coaxial  return  path  was  used  for  the  Swept  Frequency 
CW  and  Current  Pulse  tests.  The  coaxial  return  path 
provides  a  uniform  electromagnetic  field  distribution  around 
the  cylinder.  The  coaxial  return  path  is  approximately  14 
meters  long  with  a  diameter  of  2.3  meters.  The  return  path 
over  the  center  part  of  the  cylinder,  where  the  aperture  is 
located,  is  fabricated  out  of  sheet  aluminum,  2.4  meters 
long.  The  remainder  of  the  cylindrical  section  of  the 
return  path  is  constructed  out  of  wire  mesh  with  the  wires 
spaced  about  10  cm  apart  along  the  axis  and  about  5  cm  apart 
around  the  circumference  of  the  cylinder.  The  wire  spacing 
for  the  two  conical  sections  of  the  coaxial  return  path  is 
about  20  cm  along  the  cylinder  axis.  Figure  7  illustrates 
the  Lightning  Test  Cylinder  with  the  coaxial  return  path.  A 
flat  ground  plate  return  path,  however,  was  used  for  the 
Shock -Exc i tat i on  test.  The  flat  ground  plate  return,  rather 
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re  7.  Lightning  Test  Cylinder  with  Coax 


than  the  coaxial  return,  was  used  to  avoid  possible  arcing 
between  the  cylinder  and  return  path  by  the  high  voltage 
Marx  generator  during  the  Shock-Excitation  test. 

Sensors  and  Data  Acquisition 

The  electric  field  (D-dot)  and  magnetic  field  (B-dot) 
sensors  used  in  the  lightning  simulation  tests  are 
derivative  sensors,  which  respond  to  the  time  rate-of -change 
of  the  corresponding  field.  The  field  sensors,  manufactured 
by  EG&G,  are  passive  devices,  requiring  no  external  power. 
The  pertinent  equation  characterizing  the  response  of  the 
B-dot  sensors  is 

V0  =  AEQ (dB/dt)  (12) 

For  the  D-dot  sensor,  the  response  is  given  as 

V0  =  R  Aeq (dD/dt)  (13) 


where 

vo 

AEQ 

R 

B 

D 


output  voltage  (volts) 
sensor  equivalent  area  (m2) 

sensor  characteristic  load  impedance  (ohms) 
magnetic  flux  density  (webers/m2) 
electric  flux  density  (coulombs/m2) 
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In  addition  to  the  electric  and  magnetic  field  measurements. 


the  current  levels  on  the  cylinder  are  also  measured  by  a 
current  transformer  for  the  Current  Pulse  test  and  by  a 
current  shunt  for  the  Shock-Excitation  test.  The 
specifications  for  the  sensors  used  in  this  thesis  are 
provided  in  Appendix  B. 

The  output  response  of  each  sensor  is  fed  into  a  fiber 
optic  transmitter,  which  convert  the  sensor's  electrical 
response  to  an  optical  signal.  The  fiber  optic  transmitters 
are  battery  operated  and  are  activated  pneumatically . 

Because  of  the  severe  electromagnetic  environment  generated 
during  the  lightning  simulation  tests,  the  fiber  optic  links 
were  used  to  eliminate  common  mode  interference  and  ground 
loop  problems  (2:107).  A  total  of  four  fiber  optic  links 
were  used.  The  fiber  optic  transmitters  themselves  were 
shielded  to  reduce  noise  pickup.  Noise  checks  were 
performed  for  each  simulation  setup  to  insure  that  the 
transmitters  were  not  affected  by  any  noise,  such  as 
electromagnetic  interference  from  surrounding  sources. 

Fiber  optic  cables  linked  the  transmitters  to  fiber 
optic  receivers  located  in  an  EM-shielded  instrumentation 
room.  The  receivers  converted  the  optical  signals  back  to 
electrical  sensor  responses.  The  sensor  responses  are  fed 
into  two  2-channel  Tektronix  7612D  waveform  digitizers.  The 
digitizers,  which  are  triggered  by  one  of  the  sensors,  are 
set  to  read  the  data  at  2048  samples  with  a  sampling 


interval  of  5  nanoseconds .  This  provides  a  10.24 


microsecond  data  window  with  a  frequency  response  of  100 
MHz.  The  sampling  rate  is  identical  to  the  digital  data  of 
the  airborne  lightning  measurement  program.  The  data  from 
the  digitizers  were  then  read  into  a  Digital  PDP  11/34 
computer  and  recorded  on  9-track  magnetic  tape.  The  PDP 
11/34  computer  controls  data  acquisition  and  subsequently 
performs  data  reduction.  A  block  diagram  of  the  data 
acquisition  system  is  depicted  in  Figure  8.  A  listing  of 
the  computer  program  used  for  data  acquisition  is  provided 
in  Appendix  C. 

Swept  Frequency  Continuous  Wave  Test 

The  SFCW  test  measurements  were  conducted  by  Captain 
Randy  J.  Jost  of  AFIT/ENG.  Some  of  the  results  of  his 
measurements  are  incorporated  into  this  thesis.  The  setup 
of  the  SFCW  test  is  depicted  in  Figure  9.  With  a  radius 
ratio  of  2.3,  the  cylinder  and  the  coaxial  return  path 
represent  a  coaxial  transmission  line  with  a  characteristic 
impedance  of  approximately  50  ohms.  A  50  ohm  matched  load, 
therefore,  is  placed  between  the  cylinder  and  the  return 
path  in  order  to  eliminate  reflected  waves. 

A  Hewlett-Packard  3577A  network  analyzer  is  the  heart 
of  the  SFCW  svst  m.  It  generates  a  CW  current,  swept  from 
DC  to  100  M':~,  wnxch  is  amplified  and  injected  through  the 
cylinder.  -V:  -7<sG  CML-7  B-dot  sensor  is  then  used  to 


Figure  8.  Block  Diagram  of  Data  Acquisition  System 


measure  the  magnetic  field  responses  to  the  swept  frequency 


CW  induced  currents.  The  sensor  responses  are  fed  back  into 
the  network  analyzer.  Through  fiber  optic  links,  the 
response  and  the  input  current  are  recorded  by  the  Digital 
PDP  11/34  computer  and  stored  on  magnetic  disks.  The 
locations  of  the  sensor  measurements  used  in  this  thesis  are 
depicted  in  Figure  10.  The  SFCW  tests  were  conducted  on  the 
cylinder  with  the  open  aperture,  the  solid  aluminum  panel, 
the  composite  #1  panel,  and  the  composite  #2  panel. 

Current  Pulse  Measurements 

The  second  part  of  the  experimental  measurements  were 
made  using  the  current  pulse  lightning  simulation  technique. 
The  Lightning  Test  Cylinder  was  configured  with  the  coaxial 
return  path,  which  is  identical  to  the  swept  frequency  CW 
setup  except  for  the  50  ohm  load  between  the  cylinder  and 
coaxial  return.  The  setup  for  the  current  pulse  measurement 
is  illustrated  in  Figure  11. 

The  current  pulse  is  produced  by  a  generator  consisting 
of  two  banks  of  capacitors,  each  with  a  capacitance  of  eight 
microfarads.  Each  capacitor  bank  can  be  charged  to  lOOkV 
resulting  in  a  200  kV  potential  difference  capability  for 
the  current  pulse  generator.  For  this  thesis,  three 
different  current  pulse  waveforms  were  injected  into  the 
cylinder.  The  first  type  of  current  pulse  was  a  unipolar, 
doub le -exponent  ia  1  waveform  with  a  20  kA  current  peak.  A 
typical  unipolar  waveform  that  was  used  in  the  testing  is 
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amplitude  of  100  kA  is  injected  into  the  cylinder.  Figure 
13a  illustrates  a  typical  oscillatory  waveform  for  a  100 
microsecond  window.  An  oscillatory  current  waveform  with  a 
peak  amplitude  of  20  kA,  shown  in  Figure  13b,  was  also 
injected  into  the  cylinder  for  comparison  with  results  from 
the  unipolar,  douli  -exponential  current  pulse  with  the  20 
kA  peak . 

The  current  pulse  was  measured  with  a  T&M  current 
transformer  located  at  the  rear  of  the  cylinder.  The 
current  transformer  has  a  load  of  0.005  ohms.  External  and 
internal  field  measurements  of  the  cylinder  were  also  taken. 
The  internal  sensors  used  during  the  current  injection  tests 
were  EG&G  hollow  spherical  dipole  (HSD)  D-dot  sensors  for 
electric  field  measurements  and  EG&G  cylindrical  moebius 
loop  ( CML )  B-dot  sensors  for  magnetic  field  measurements. 
External  magnetic  field  measurements  for  the  unipolar  20  kA 
current  pulse  tests  were  also  performed  by  CML  B-dot 
sensors.  External  magnetic  field  measurements  for  the 
oscillatory  current  waveforms  were  done  by  multi-gap  loop 
MGL )  B-dot  sensors.  Asymptotic  conical  dipole  ( ACD )  D-dot 
sensors  were  utilized  for  external  electric  field 
measurements.  Figure  14  depicts  the  locations  of  the  D-dot 
and  B-dot  sensors  on  the  cylinder  where  measurements  were 
recorded  during  the  current  pulse  tests.  The  current  pulse 
tests  were  conducted  with  the  aperture  open,  or  covered  by 
the  solid  aluminum,  the  composite  #1,  and  the  composite  # 2 
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Figure  13.  Oscillatory  Current  Waveforms 
(a)  100  kA  Peak  and  (b)  20  kA  Peak. 
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Sensor  Locations  for  Current  Pulse  Test 


panels  using  the  20  kA  unipolar,  20  kA  oscillatory,  and  100 
kA  oscillatory  injected  waveforms. 

Shock-Exc i tat i on  Measurements 

The  final  lightning  simulation  test  that  was  conducted 
on  the  cylinder  incorporated  the  shock -exc i tat i on  technique. 
The  test  setup  is  depicted  in  Figure  15,  where  a  flat  ground 
plate,  rather  than  a  coaxial,  return  path  was  used  in  the 
configuration.  The  cylinder  was  supported  75  centimeters 
above  the  ground  plate  by  two  dielectric  stands. 

The  excitation  source  is  provided  by  a  Marx  generator. 
The  Marx  generator  used  in  this  shock-excitation 
configuration  consists  of  29  stages  of  0.7  microfarads 
capacitors.  Each  capacitor  can  be  charged  to  50  kV  which 
results  in  a  total  discharge  voltage  potential  of  about  1.45 
megavolts.  For  the  shock -exc i tat i on  measurements  in  this 
thesis,  the  Marx  generator  was  charged  to  deliver  one 
megavolt.  The  waveshaping  network  on  the  discharge  end  of 
the  Marx  generator  consists  of  a  total  resistance  of  1800 
ohms  in  series  with  an  inductor  of  116  microhenries.  A 
grounded  electrostatic  shield  is  placed  between  the  Marx 
generator  and  the  test  cylinder  in  order  to  shield  the 
cylinder  from  the  fields  generated  by  the  Marx  as  it 
discharges  . 

The  input  and  output  spark  gaps  were  separated  at 
various  distances  to  examine  its  effects  on  charging  and 


discharging  transients  on  the  cylinder. 


Input  spark  gaps  of 


Shock-Excitation  Test  Setup. 
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Figure  25 .  P3  Exte 
(a)  20  kA  Unipolar, 
Osc i 1 latory . 


electric  field  as  the  current  is  injected  into  the  cylinder. 
However,  the  positive  field  following  the  negative  electric 
field,  which  is  seen  in  the  unipolar  results,  is  not 
apparent  in  the  external  responses  of  the  20  kA  and  100  kA 
oscillatory  currents.  The  electric  field  responses 
oscillate  in  the  negative  region  for  almost  the  entire  10.24 
microsecond  window. 

The  external  magnetic  field  responses  (Figure  25)  tend 
to  follow  tne  same  waveshape  as  the  corresponding  injected 
current  pulse  --  unipolar  and  oscillatory.  This  is  due  to 
the  direct  proportional  relationship  from  Ampere's  law 
between  the  H-field  and  current  for  a  cylindrical  conductor. 
Although  the  injected  excitation  source  is  actually  a 
negative  current  pulse,  the  H-field  responses  are  plotted 
for  convenience  to  indicate  a  positive  current. 

External  Field  Responses  to  Shock-Excitation 

The  shock -exc i tat i on  responses  presented  in  this 
section  is  for  solid  aluminum  panel  configuration  using  the 
60  centimeter  input  spark  gap  and  the  30  centimeter  output 
spark  gap  setup.  Figure  26  shows  the  external  electric 
field  response  for  the  s hock -e xc i ta t i on  simulation.  Note 
the  initial,  almost  linear,  negative  increase  in  the 
electric  field  due  to  input  spark  gap  breakdown,  which  is 
followed  by  charging  phase  of  the  simulation.  Once  the  arc 


across  the  output  spark  gap  has  been  established,  the 


IE  3  ( V/M )  E(A1 ) 


Figure  26.  P3  External  E-Field  Response  for 
Shock -Excitation . 


electric  field  immediately  oscillates  with  the  excitation 
current  flowing  through  the  cylinder. 

The  external  H-field  response  is  illustrated  in  Figure 
27.  During  the  input  spark  gap  breakdown  and  the  charging 
phase  of  the  shock -exc i tat  ion  simulation,  the  H-field  has  a 
negative  rate  of  change.  The  H-field  then  experiences  a 
very  rapid  positive  change  once  the  output  arc  has  been 
established.  Similar  to  the  E-field  response,  the  H-field 
then  follows  the  excitation  source  as  the  oscillating 
current  is  discharged  through  the  cylinder. 

Comparison  cf  Transfer  Functions 

This  section  will  discuss  the  comparisons  between  the 


transfer  functions  that  were  derived  from  the  SFCW,  current 
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Figure  27.  P3  External  H-Field  Response  for 
Shock -Excitation. 


pulse,  and  shock-excitation  lightning  simulation  techniques. 
Figure  28  depict  the  transfer  function  of  the  external 
H-field  determined  by  the  SFCW  technique  for  the  P3  setup. 

A  corresponding  linear  plot  of  the  transfer  function  is  also 
shown  in  order  to  clearly  identify  particular  frequencies. 
The  SFCW  derived  transfer  functions  are  first  compared  to 
those  obtained  from  the  current  pulse  tests,  which  are  shown 
in  Figures  29  through  31,  for  the  external  H-field 
responses.  In  the  case  of  the  current  pulse  technique,  one 
notes  the  decreasing  magnitude  of  the  transfer  functions  of 
the  external  magnetic  field  response  to  the  excitation 
current  at  higher  frequencies.  This  general  trend  is  seen 
for  both  the  unipolar  current  pulse  (Figure  29)  and 
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Figure  28  . 
for  SFCW . 


P3  External  H-Field  Transfer  Function 


Figure  29.  P3  External  H-Field  Transfer  Function 
for  20  kA  Unipolar  Current  Pulse. 


oscillatory  current  pulse  results  (Figures  30  and  31).  This 
is  unlike  the  transfer  functions  obtained  from  the  SFCW 
technique  which  have  a  relatively  constant  magnitude  until 
about  80  MHz.  The  oscillatory  transfer  functions,  20  kA  and 
100  kA,  have  a  faster  decrease  in  magnitude  with  increasing 
frequency,  however,  when  compared  to  the  unipolar  case.  The 
differences  between  the  transfer  functions  of  the  various 
simulation  techniques  are  also  seen  when  they  are  compared 
to  those  of  the  shock-excitation  technique.  Figure  32  shows 
the  external  magnetic  field  transfer  function  that  is 
derived  with  respect  to  the  excitation  current.  For  the 
shock -exc i tat i on  case,  the  transfer  function  tend  to 
increase  in  magnitude  at  higher  frequencies.  This  effect  is 
possibly  due  to  the  very  fast  risetimes  of  the  excitation 
current  that  are  obtained  with  the  shock-excitation 
technique . 

The  linear  plots  of  the  transfer  functions  clearly 
identify  specific  lobes.  In  Figure  28,  spikes  in  the 
transfer  function  occur  at  17,  29,  41,  52,  64,  and  87  MHz 
for  the  SFCW  case.  The  most  prominent  lobe  appears  at  41 
MHz  for  the  solid  aluminum  panel  conf igurat ion  P3. 

The  external  H-field  transfer  functions  for  the  current 
pulse  case  (Figures  29  through  31)  indicate  lobes  at  6,  8, 
17,  and  20  MHz,  in  addition  to  a  number  of  lobes  between  40 
to  60  MHz,  for  both  the  unipolar  and  the  oscillatory  current 
pulses.  Those  lobes  occuring  at  17  MHz  and  between  the  40 
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to  60  MHz  bandwidth  can  be  correlated  to  lobes  at  the  same 
frequencies  for  the  SFCW  derived  transfer  function. 

The  external  H-field  transfer  function  for  the 
shock-excitation  setup  (Figure  32)  shows  very  large  spikes 
at  57  and  67  MHz.  These  character istic  lobes  are  possibly 
due  to  the  facility  effects  such  as  the  return  path  and  the 
generator  used  during  the  shock-excitation  tests. 

The  external  E-field  transfer  functions  resulting  from 
the  current  pulse  measurements  are  illustrated  in  Figures  33 
through  35.  These  transfer  functions  show  a  relatively 
constant  magnitude  until  the  first  large  lobe  occuring  at  6 
MHz.  All  of  the  transfer  functions  then  decrease  in 
magnitude  for  frequencies  beyond  6  MHz.  Unlike  the  transfer 
function  of  the  20  kA  unipolar  (Figure  33)  and  the  100  kA 
oscillatory  (Figure  35)  current  pulses,  the  transfer 
function  of  the  20  kA  oscillatory  pulse  (Figure  34) 
experiences  a  number  of  nulls  prior  to  the  6  MHz  mark. 
However,  in  all  cases,  unipolar  and  oscillatory,  another 
lobe  is  present  at  8  MHz  that  is  soon  followed  by  nulls  at 
10  and  16  MHz.  Other  lobes  occur  at  13,  17,  and  20  MHz. 

The  shock-exc i tat i on  transfer  function  of  the  external 
electric  field  is  depicted  in  Figure  36  for  the  solid 
aluminum  panel  conf igurat i on .  As  the  frequency  increases, 
the  magnitude  of  the  transfer  functions  also  increase  until 
the  first  relative  null  at  approximately  600  kHz.  The 
magnitude  then  decreases  as  the  frequency  approaches  12  MHz. 


Figure  34.  P3  External  E-Field  Transfer  Function 
for  20  kA  Oscillatory  Current  Pulse. 


The  shock-excitation  transfer  function  of  the  external 


E-field  has  character ist ic  lobes  at  400  and  800  kHz  and  a 
null  at  12  MHz.  At  the  higher  frequencies,  prominent  spikes 
appear  at  57,  69,  and  72  MHz,  as  well  as  the  extremely  large 
lobe  at  67  MHz.  Those  spikes  appearing  at  57  and  67  MHz 
also  occur  in  the  shock-excitation  transfer  function  for  the 
external  H-field  (Figure  32).  Again,  these  two  spikes  may 
be  an  effect  resulting  from  the  facility  setup  of  the 
shock-excitation  simulation;  spec i f ical ly,  the  flat  ground 
plate  return  path  and  the  Marx  generator. 

Threat-Level  Field  Responses 

Since  the  transfer  functions  that  were  derived  for  each 
of  the  various  simulation  techniques  uniquely  describe  the 
induced  electric  and  magnetic  fields  on  the  test  cylinder, 
the  threat-level  field  responses  can  be  predicted.  Using 
the  assumption  that  the  electromagnetic  interaction  is  a 
linear  process,  the  threat-level  field  responses  are 
determined  by  folding  the  threat-level  current  waveform  into 
the  transfer  function.  The  threat  waveform  used  in  this 
thesis  for  response  prediction  is  the  Current  Waveform 
Component  A,  which  represents  the  severe  case  of  a  first 
return  stroke.  This  double  exponential  waveform  has  a  peak 
amplitude  of  200  kA  and  an  action  integral  of  2  x  102  A2 s . 


The  peak  response  values  of  the  predicted  full 
threat-level  electric  field  and  magnetic  field  are  shown  in 
Table  1.  Both  external  (Location  A1 )  and  internal  (Location 
A3)  field  responses  were  computed  for  the  various  panel 
configurations  and  lightning  simulation  techniques.  Since 
the  transfer  functions  derived  from  the  SFCW  technique  were 
determined  only  for  the  external  H-field,  the  SFCW 
threat-level  predictions  do  not  include  the  internal  and 
external  E-field  responses  and  the  internal  H-field 
responses  . 

In  comparing  the  full  threat-level  responses,  the 
results  of  the  solid  panel  configuration  are  used.  The  peak 
external  H-field  responses  for  the  different  simulation 
techniques  are  graphically  illustrated  in  Figure  38.  The 
SFCW  technique  has  the  largest  threat-level  response  with  a 
predicted  peak  value  of  58.3  kA/m.  This  is  about  twice  the 
peak  H-field  response  values  predicted  by  the 
shock -exc i tat i on  technique  and  the  20  kA  unipolar  current 
pulse  technique.  The  two  smallest  predicted  values  are 
given  by  both  oscillatory  current  pulse  methods. 

Figure  39  is  a  graphical  presentation  of  the  peak 
external  E-field  responses  to  the  threat  current  waveform. 
The  magnitude  of  the  peak  electric  field  response  predicted 
by  the  shock -exc i tat i on  technique  is  over  100  times  the  peak 
response  values  predicted  by  the  unipolar  and  oscillatory 
current  pulse  techniques.  The  shock -exc  i  tat i on  E-field  peak 


Table  1 


Responses  to  Threat-Level  Current  Waveform 
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Field 

PO 

PI 

P  2 

P  3 

SFCW 

H  (  A1 ) 

l  A/ml 

31 . 2E3 

52.6E3 

62 . 1E3 

58 . 3E3 

C20u 

H(  A1  ] 
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H  .  1E3 

26 . 6E3 

27 . 2E3 

26 . 0E3 
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(A/ml 

5  .  7E3 

17 . 3E3 

20 . 8E3 

10 . 2E3 

ClOOo 

H(A1) 

(  A/ml 

3 . 7E3 

10 . 8E3 

12  .  6E3 

11 . 4E3 

SE 

H  (  A1  ) 

(  A/m] 
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496E3 

431E3 

537E3 

48SE3 

C20o 

E(A1) 

(  V/mi 

86E3 

37E3 

8  7E3 

1 1  7  E  3 

ClOOo 

E  (  A1  ) 

[  V/m  1 

76E3 
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52 . 9E6 
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H  (  A3  ) 
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6  .  6E3 
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2  .  3 

C20o 

H  (  A3  ) 

1  A/m) 

3.  2E3 

432 

217 

1  .  1 

ClOOo 

H  (  A3  ) 

(  A/m  1 

2  .  3E3 

400 

226 

0  .  4 

SE 

H(  A3  I 

1  A/m  1 

21 . 1E3 

916 

321 

31.1 

C20u 

E(  A3  ) 

(  V/m  1 

4  4  3  E  3 

1  4  E  3 

6  E  3 

3  E  3 

C20o 

E(  A3  ) 
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6  1 E  3 

2 1 E  3 

1  4E3 

3E3 

ClOOo 

E(  A3  ) 

(  V/m  1 

7  4  E  3 

1  3  E  3 

9E3 

1  E  3 

'  E 

E  (  A3  ) 

( V/ml 

26 . 9E6 

7  1  E  3 

4  3E3 

5  7  E  3 

Notation: 

3  FCW 

=* 

swept  frequency  continuous 

w  a  ve 

C  2  0  o 

= 

20  kA  unipolar  current  pulse 

C  2  Oo 

= 

20  kA  oscillatory  current 

pulse 

ClOOo 

= 

100  kA  oscillatory  current 

pulse 

SE 

= 

shock-excitation 

P0 

= 

open  aperture 

FI 

= 

composite  11  panel 

P  2 

= 

composite  12  panel 

P  3 

= 

solid  aluminum  panel 

is  49.7  MV/m;  whereas,  the  20  kA  unipolar  current  pulse 
response  peak  has  a  magnitude  of  485  kV/m.  As  in  the 
external  H-field  threat-level  responses,  the  predictions  by 
the  two  oscillatory  current  pulse  methods  have  the  smallest 
magnitudes  of  all  of  the  test  techniques.  Both  the  20  kA 
and  the  100  kA  oscillatory  current  pulse  transfer  functions 
compute  an  external  E-field  peak  of  117  kV/m  from  the  200  kA 
threat  current  waveform. 

The  external  H-field  and  E-field  predictions  are  also 
extended  to  the  internal  field  responses.  The  peak  internal 
H-field  and  E-field  responses  to  the  threat  current  waveform 
are  illustrated,  respectively,  in  Figure  40  and  Figure  41 
for  the  solid  aluminum  panel  configuration.  In  both  cases, 
the  peak  internal  field  responses  are  dominated  by  the 
prediction  generated  by  the  shock-excitation  simulation 
technique.  The  20  kA  unipolar  case  produces  the  largest 
internal  H-field  and  E-field  responses  of  the  three  current 
pulse  techniques. 

The  predicted  full  threat-level  responses  can  also  be 
compared  to  the  actual  measured  responses  of  the  current 
pulse  and  the  shock -exc i tat i on  lightning  simulation 
techniques.  The  measured  peak  magnetic  and  electric  field 
responses  for  the  various  panel  configurations  are  shown  in 
Table  2.  When  compared  to  the  predicted  full  threat-level 
responses  shown  in  Table  1,  one  immediately  observes  that 
the  actual  measured  peak  field  responses  for  the  100  kA 


psfvrv; 


oscillatory  current  pulse  tests  are  very  close  to  the 
corresponding  predicted  threat-level  values.  This 
observation  appears  to  be  true  for  all  four  panel 
configurations  for  both  external  and  internal  field 


responses 


magnetic  and  electric. 


In  this  section,  the  measured  field  responses  and  the 

transfer  functions  resulting  from  the  current  pulse  and  the 

shock-excitation  simulation  techniques  are  compared  to 

measured  lightning  data  from  the  lightning  characterization 

program  during  1984  and  1985.  The  maximum  levels  of 

electric  and  magnetic  flux  densities  measured  by  the  CV-580 

aircraft  during  the  two  year  program  were  22  Coul/m^s  and 

3950  Wb/m2/s,  respectively  {18: iv).  Assuming  that  the 

relative  permittivity  and  relative  permeability  are  equal  to 

12 

one,  this  would  correspond  to  an  E-dot  (dE/dt)  of  3.7  x  10 
V/m/s  and  an  H-dot  (dH/dt)  of  3.1  x  10®  A/m/s. 

The  peak  E-dot  and  H-dot  responses  for  the  current 
pulse  and  shock-excitation  techniques  are  shown  in  Table  3 
for  the  solid  aluminum  panel  configuration.  As  can  be  seen 
from  this  table,  both  the  dE/dt's  and  dH/dt's  generated  by 
the  current  pulse  and  the  shock-excitation  techniques  are 
very  close  to  those  measured  in  the  airborne  case.  The  peak 
H-dot  responses  achieved  by  both  the  current  pulse  and 
shock -exc i tat i on  methods  appear  to  be  sufficient  in 
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simulating  this  airborne  event  parameter.  However,  it 


Table  3 


ISiSmSrCT 


Test 

Peak 

H-Dot(Al)  [A/ m/s] 

Peak 

E-Dot (Al)  [V/m/s] 

C20u 

10.6  x  109 

1.2  x  1012 

C20o 

3.8  x  109 

0.7  x  1012 

ClOOo 

12.6  x  10® 

1.2  x  1012 

SE 

9.0  x  10® 

4.5  x  1012 

In-Flight 

3.1  x  10® 

3.7  x  1012 

Notation : 
C20u 
C20o 
ClOOo 
SE 


20  kA  unipolar  current  pulse 
20  kA  oscillatory  current  pulse 
100  kA  oscillatory  current  pulse 
s hock -excitation 
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Figure  42.  Transfer  Function  of  External  H-Field 
Response  for  CV-580  Airborne  Lightning  Measurement 


16  MHz  and  reflects  a  full  wavelength  dimension  of  the 
coaxial  return  path.  The  test  cylinder  itself  is 
represented  similarly  by  the  lobes  occuring  at  6  and  20  MHz 
which  correspond  approximately  to  40  and  12  meters, 
respectively.  The  two  prominent  spikes  for  the 
shock-excitation  case,  shown  in  Figure  32,  represent 
appr oximate ly  4  meters  for  the  spike  at  57  MHz  and  3.6 
meters  for  the  spike  at  67  MHz. 
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Conclusions  and  Recommendations 


Based  on  the  experimental  results  obtained  during  this 
thesis  and  presented  in  Chapter  4,  several  conclusions  can 
be  drawn  which  answers  the  questions  initially  posed  for 
this  research  effort.  One  of  the  goals  of  this  thesis  was 
to  compare  the  transients  produced  by  the  current  pulse  and 
shock-excitation  simulation  methods  with  swept  frequency 
continuous  wave  measurements.  The  comparisons  were  made 
with  respect  to  external  H-field  transients  that  were 
linearly  scaled  to  the  full  threat-level  current  waveform. 
The  results,  shown  in  Figure  38,  indicate  that  the  SFCW 
linear  scaled  measurement  has  the  largest  predicted  H-field 
transient  of  all  of  the  test  -techniques  investigated  in  this 
thesis.  The  results  from  both  oscillatory  current  pulse 
methods  have  the  lowest  predicted  H-field  responses.  From 
these  results,  one  can  see  the  problems  of  linearly  scaling 
low-level  transients.  The  SFCW  method,  which  uses  very  low 
energy  levels  of  approximately  25  watts,  appear  to 
overestimate  the  full  threat-level  response.  The  problem 
seems  to  be  the  extreme  extrapolation  that  is  necessary  to 
scale  the  SFCW  measurement  to  energy  levels  that  are 
associated  with  a  threat  current  waveform  with  a  peak 
amplitude  of  200  kA.  This  problem  can  also  be  attributed  to 
the  very  large  predicted  H-field  transient  for  the 
shock -exc i tat i on  method.  The  measured  peak  current  of  the 


shock-excitation  method  is  approximately  1.4  kA,  which  is 
then  linearly  scaled  to  200  kA.  On  the  other  hand,  the 
external  H-field  measurement  from  the  100  kA  oscillatory 
current  pulse  method  requires  a  relatively  small 
extrapolation  to  the  full  threat-level  current  waveform;  and 
therefore,  it  appears  to  predict  more  reasonable  results. 
This  can  be  seen  from  the  actual  measured  transients  for  the 
100  kA  oscillatory  current  pulse  case  in  Table  2.  From 
Table  2,  one  sees  that  the  actual  measured  results  for  the 
100  kA  oscillatory  current  pulse  are  very  close  to  the 
predicted  full  threat-level  values  due  to  the  smaller 
extrapolation  that  is  required.  These  results  point  out  the 
validity  of  the  assumption  of  a  linear  system  because  it 
tends  to  predict  conservatively  higher  transient  values. 
However,  it  also  points  out  that  extreme  care  must  be 
exercised  when  analyzing  these  linearly  extrapolated  values, 
since  scaling  the  measured  responses  overestimates  the 
predicted  values.  The  less  the  extrapolation  that  is 
necessary,  the  more  accurate  the  predicted  induced  transient 
values  will  be  ( 5 ) . 

In  analyzing  the  experimental  time  rate-of -change 
results  of  the  current  pulse  lightning  simulation  tests,  one 
sees  that  the  dH/dt  responses  produced  by  the  current  pulse 
tests  are  sufficient  in  replicating  the  airborne  measurement 
of  3.1  x  10*  A/m/s.  However,  the  current  pulse  method  was 
not  able  to  achieve  the  dE/dt  responses  of  the  airborne 


case,  which  is  3.7  x  lO12  V/m/s .  The  difference  is  small; 
the  dE/dt  transients  are  below  the  airborne  event  by 
approximately  a  factor  of  3.  It  appears  that  only  the 
shock-excitation  technique  was  able  to  generate  the  peak 
dE/dt  responses  that  is  sufficiently  high  enough  to  exceed 
the  CV- 580  airborne  lightning  strike  measurement. 

In  terms  of  electric  field  transients,  the  in-flight 
measurements  of  the  CV-580  recorded  maximum  E-field 
transients  of  about  200  kV/m  (18:154).  From  results  of  the 
peak  external  E-field  responses  shown  in  Table  2,  only  the 
shock -exc i tat i on  technique  is  able  to  simulate  airborne 
E-field  transient.  The  shock-excitation  test  method 
produces  a  peak  E-field  transient  of  241  kV/m  for  the  solid 
aluminum  panel  configuration  which  exceeds  the  measured 
airborne  response. 

Therefore,  in  order  to  generate  the  electric  field  and 
the  time  rate -o f -change  of  the  electric  field  that  are 
measured  for  the  aircraft/lightning  interaction  event,  it 
appears  to  be  necessary  to  employ  the  s hock -exc i tat i on  test 
method.  However,  these  values  for  the  electric  field  and 
the  time  rate-of -change  of  the  electric  field  can  also  be 
achieved  through  the  current  pulse  simulation  technique  by 
linearly  extrapolating  the  measured  responses.  The  full 
threat-level  responses,  as  shown  in  Table  1,  are  obtained  in 
this  thesis  by  deriving  the  transfer  function  of  the  field 
response  with  respect  to  the  injected  excitation  current. 


and  then  folding  the  full  threat-level  current  waveform  into 
the  transfer  function.  It  is  important  that  one  takes  into 
account  the  degree  of  over est i mat i on  that  occurs  when 
linearly  scaling  the  responses  to  the  full  threat-level. 

The  amount  by  which  the  induced  transients  are  to  be 
extrapolated  in  the  analysis  should  be  minimized.  One 
option,  which  may  be  impractical  to  use  on  large  aircraft, 
is  to  incorporate  very  high  current  generators  in  the 
current  pulse  simulation  tests. 

The  effects  that  the  current  pulse  test  method  is  not 
able  to  simulate  are  the  effects  of  nearby  lightning  and  the 
effects  of  the  approaching  stepped  leader  of  a  lightning 
strike.  However,  the  shock-excitation  technique  simulates 
these  effects,  respectively.  Cry  the  rapidly  changing  E-field 
at  the  input  spark  gap  and  by  the  charging  phase  of  the 
simulation.  The  electric  field  responses  generated  by  the 
shock-excitation  technique  are  very  similar  to  those 
measured  by  the  CV-580  in-flight  tests.  A  typical  analog 
measurement  of  the  E-field  during  a  lightning  strike  to  an 
aircraft,  which  is  shown  in  Figure  43,  resembles  very 
closely  to  the  results  of  the  charging  and  discharging 
phases  of  the  shock -exc i tat i on  tests. 

In  terms  of  current,  magnetic  field,  and  electric  field 
transients,  the  experimental  results  of  this  thesis  indicate 
that  the  current  pulse  technique  is  able  to  reproduce  the 
effects  measured  during  in-flight  lightning  strikes  to  the 


Figure  43.  Typical  E-Field  Waveform  During  a 
Lightning  Strike  to  Airborne  CV-580  Aircraft. 


CV-580  aircraft.  An  important  note  is  that  the  comparisons 
conducted  in  this  thesis  were  of  a  limited  number  of 
parameters  which  only  inc luded^the i r  peak  magnitudes  and 
their  peak  time  rate-of -changes . 

The  experimental  results  of  this  thesis,  however,  are 
not  able  to  answer  the  question  of  whether  the  the  current 
pulse  technique  alone  is  sufficient  to  assess  the 
suscept i bl i ty  and  vulnerability  of  an  aircraft  to  the 
lightning  threat.  The  question  in  this  case  is  whether  the 
prereturn  stroke  phase  of  the  lightning  strike  presents  a 
unique  hazard  to  aircraft  (4:6).  If  it  does,  then  the  next 
question  is  whether  the  prereturn  stroke  phase  can  be 
neglected  if  the  aircraft  is  tested  and  protected  using  the 


results  from  current  pulse  testing.  These  questions  on 
lightning  simulation  testing  require  further  investigation. 


Should  the  current  pulse  method  prove  to  be  inadequate,  then 
the  shock-excitation  technique  can  be  incorporated  to 
complement  lightning  simulation  testing. 

Re commend at  ions 

This  thesis  effort  took  a  first  step  in  assessing  the 
capabilities  of  lightning  simulation  techniques  by  comparing 
the  electrical  parameters  to  airborne  lightning  strike 
measurements.  There  are  several  areas  that  should  be 
investigated  in  future  studies  lightning  simulation. 

Further  comparison  and  analysis  of  the  external  and  internal 
measurements  of  the  electric  and  magnetic  field  responses 
should  be  examined  for  the  various  lightning  simulation 
techniques.  This  study  can  include  analysis  of  the  electric 
and  magnetic  field  responses  for  different  composite 
aircraft  structures. 

A  possible  follow  on  thesis  in  assessing  lightning 
simulation  would  be  to  calibrate  the  induced  transients 
produced  by  lightning  simulation  in  terms  of  the  measured 
airborne  lightning  strike  transients.  The  goal  would  then 
be  determining  a  scaling  function  that  would  correlate  the 
simulated  lightning  transients  to  the  actual  airborne  event. 
The  results  from  this  study  may  improve  the  capability  of 
lightning  simulation  technology  to  faithfully  reproduce  the 
effects  of  a  lightning  strike  to  flying  aircraft. 


This  appendix  discusses  the  characterization  program 


where  transient  responses  were  recorded  on  a  specially 
instrumented  lightning  research  aircraft  during  an  actual 
lightning  strike.  During  1984  and  1985,  AFWAL/FI ESL , 
working  with  the  FAA,  U.S.  Navy,  NASA,  and  ONERA,  conducted 
an  airborne  lightning  characterization  program.  The  program 
measured  and  recorded  skin  current  distributions  and 
electromagnetic  fields  resulting  from  direct  lightning 
attachments  to  a  flying  aircraft.  A  CV-580  aircraft, 
provided  by  the  FAA,  was  instrumented  by  AFWAL/FIESL  and 
flown  in  and  near  active  Florida  thunderstorms  at  altitudes 
between  1,800  and  18,000  feet-. 

Aircraft  Instrumentation 

The  CV-580  was  instrumented  with  a  total  of  27  sensors 
(16:2).  The  instrumentation  of  interest  for  this  thesis  are 
the  externally  mounted  current  shunts,  electric  field 
sensors,  and  magnetic  field  sensors.  The  current  shunts 
were  mounted  on  the  right  wing  tip  (IRW)  and  the  left  wing 
tip  (ILW)  of  the  CV-580  aircraft  (1:2).  For  the  test 
flights  during  1985,  additional  current  shunts  were  mounted 
on  top  of  the  vertical  stabilizer  (IVS)  and  in  a  tail  boom 
(ITB)  behind  the  vertical  stabilizer  (1:2). 

Derivative  magnetic  field  (surface  current)  sensors 
were  mounted  on  top  of  the  forward  fuselage  (JSFF),  the  aft 
fuselage  (JSAF),  and  the  left  wing  between  the  fuselage  and 


the  engine  pod  (JSTLW).  Surface  current  sensors  were  also 
located  between  the  fuselage  and  the  engine  pods  on  the 
bottom  of  the  left  (JSLW)  and  right  (JSRW)  wings. 

Derivative  electric  field  (displacement  current)  sensors 
were  mounted  on  top  of  the  forward  fuselage  (JNFF)  and  on 
top  of  the  right  wing  (JNTRW),  between  the  fuselage  and  the 
engine  pod.  Derivative  electric  field  sensors  were  also 
placed  on  the  bottom  of  the  right  (JNRW)  and  left  (JNLW) 
wing  tips  (1:2).  Figure  44  shows  the  placement  of  the 
various  sensors  on  the  CV-580  aircraft. 

The  high  frequency  data  from  the  sensors  were  recorded 
by  six  Textronix  7612D  waveform  digitizers,  each  with  two 
input  channels.  The  digitizers  were  set  to  record  2048 
samples  at  5  nanosecond  intervals.  This  provided  a 
recording  window  of  10.24  mi-er oseconds  of  digital  data  with 
an  upper  bandwidth  of  100  MHz.  The  low  frequency  data,  DC 
to  500  kHz  in  the  FM  channels  and  400  Hz  to  2  MHz  in  the 
direct  channels,  were  recorded  continuously  by  a  28-channel 
Honeywell  101  analog  recorder. 

Ground  Station 

Ground  station  measurements  were  also  conducted  during 
the  lightning  character i zat i on  program.  A  ground  station 
trailer,  which  was  located  north  of  Kennedy  Space  Center, 
was  instrumented  with  two  flush  plate  electric  field 
antennas  and  two  crossed-loop  magnetic  field  sensors  (16:2' 


As  presented  in  Chapter  3,  various  magnetic  and 
electric  field  sensors  were  used  during  the  lightning 
simulation  measurements  of  this  thesis.  This  appendix 
contains  the  data  sheets  of  the  specifications  for  the  EG&G 
sensors  used.  The  sensor  data  sheets  include  the 
Cylindrical  Moebius  Loop  (C ML),  Multi-Gap  Loop  ( MGL ) ,  Hollow 
Spherical  Dipole  ( HSD ) ,  Flush  Plate  Dipole  ( FPD ) ,  and 
Asymptotic  Conical  Dipole  (ACD)  sensors. 


SPECIFICATIONS 


Parameter 

CML-7  or  CML-S7 

A  (m2) 

2  x  10'2 

eq 

Frequency  Response  (3  dB  point) 

38  MHz 

Risetime  (10-90%) 

9  ns 

Maximum  Output 

4  kV  (with  TCC) 

Maximum  Field  Change 

2  x  10®  tesla/sec 

Output  Connector 

TCC* 

Mass 

1.  0  kg 

Dimensions  (cm) 

L 

35.6 

D 

11.9 

H 

12.4 

A 

6.4 

*  100-ohm  twinaxial  connector  (Data  Sheet  1340);  Two  50-ohm  SMA 
connectors  optional. 


“  5 


OUTPUT  I 

CONNICTO*  , 


(Data  and  Specifications  Subject  to  Change  without  Notice) 
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SPECIFICATIONS 


Parameter 

HSD-2 (R) 

HSD-4  (R) 

A  (m2) 
eq 

Frequency  Response 

1  x  10'1 

1  x  10'2 

>45  MHz 

>  150  MHz 

(3  dB  Point) 

Risetime  (Tp  10_90) 

s  7.  4  ns 

s  2.  3  ns 

Maximum  Output 

5  kV 

5  kV 

Output  Connector 

TCC  -  100  ohm* 

TCC  -  100  ohm* 

Mass 

2.15kg 

0.  59  kg 

Dimensions  (cm) 

L 

62.  5 

44.7 

W 

28.  5 

8.  9*» 

H 

20.6 

6.9 

t 

0.3 

0.16 

♦100-ohm  Twinaxial  Connector  (EG&G  Data  Sheet  1340) 
♦♦Connector  width  is  greater  than  W  on  Model  4  (11. 7  cm) 


NOTES: 


During  use,  this  sensor  must  be  supported  by  dielectric 
materials  and  positioned  at  least  two  sensor  diameters 
from  any  conducting  surfaces. 

A  DLT-96  balun  (EG&G  Data  Sheet  1300)  can  be  used  to 
transform  100-ohm  balanced  output  to  50-ohm  unbalanced 
output  for  telemetry  and  recording. 


(Data  and  Specifications  Subject  to  Change  without  Notice) 


SPECIFICATIONS: 


FPD-1 


FPD-2 


Aeq  <m2) 

1  x  10‘2 

2  x  10'2 

Frequency  Response 
(3  dB  point) 

>  350  MHz 

>70  MHz 

Risetime  (T.  1(J_9(J) 

c  1  ns 

<  5  ns 

Maximum  Output 

5  kV 

5  kV 

Output  Connector 

GR874L  (50  ohms) 

GR874L  (50  ohms) 

Mass 

4  kg 

2  kg 

Dimensions  (cm) 

W 

43.2 

28.  3 

t 

1.0 

0.  5 

D 

5.4 

5.2 

(Data  and  Specifications  Subject  to  Change  without  Notice) 
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SPECIFICATIONS 


ACD-Sl(R)* 

A  CD-  1(A)* 

ACD-3IA) 

ACD-S3IR1 

ACD-5IR) 

ACD-6(A) 

ACD-9(A) 

Aeq  (m  ) 

1  x  10“4 

1  x  10“* 

1  x  10~2 

1  x  10‘2 

1  x  10° 

1  x  10~’ 

1  x  10*3 

Frequency  Response 

>11.0  GHz 

>11.  0  GHz 

>1.  10  GHz 

>1.  10  GHz 

>110  MHZ 

>350  MHZ 

>3.  5  GHz 

(3-dB  point) 

Risetime  (Tr  20-90) 

<0.032  ns 

<0.032  ns 

<0.32  ns 

<0.32  ns 

<3.  2  ns 

<  1.  0  ns 

<  0.  10  ns 

Maximum  Output 

125  V 

125  V 

1000  V 

5  kV 

5  kV 

5  kV 

1000  V 

Output  Connector 

OSSM  274 

ARMM 

ARM 

GR** 

GR 

GR 

ARM 

son 

4052- 

2052- 

874L 

874L 

874L 

2052- 

coaxial 

0000 

0000 

500 

son 

50  n 

0000 

Mass 

12  g 

ISg 

327  g 

550  g 

60  kg 

1.3  kg 

60  g 

Dimensions  (cm) 

L 

7.62 

— 

— 

22.  2 

152.40 

— 

... 

W 

2.  54 

2.  54 

13.97 

14.3 

121.92 

28.2 

5.  59 

H 

0.59 

0.59 

5.18 

5.18 

49.14 

15.56 

1.79 

t 

0.  10 

0.  10 

0.  254 

0.  24 

0.84 

0.  32 

0.  236 

T 

1.32 

1.18 

1. 82 

3. 19 

4.13 

6.03 

1.82 

•  Axial  or  Radial  Outjwt  Specified  by  Designations  ACD-N(A)  or  ACD-N(R),  Respectively,  where  N  »  1  or  3. 
•*OSM  Connector  Available  on  Request. 


AXIAL 


RADIAL 


CONNECTOR  GUARO 
(ACD-6  Only) 


*  Fltiiblt  Ground  Tlint  Vinion  I 
and  Coaiial  fflra  lit  tRIa  Stdioi 


(Data  and  Specifications  Subject  to  Change  without  Notice' 
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This  appendix  contains  a  complete  listing  of  the 
computer  program,  written  by  AFWAL/FIESL,  that  was  used  for 
data  acquisition  during  the  current  pulse  and 

shock-excitation  lightning  simulation  tests  on  the  Lightning 
Test  Cylinder.  The  program  controls  the  two  2-channel 
Tektronix  7612D  waveform  digitizers  and  performs  data 
reduction.  The  data  is  stored  on  9-track  magnetic  tape. 

The  program  runs  on  the  Digital  PDP  11/34  using  the  TEK  SPS 
BASIC  Signal  Processing  Package. 


1  REM  DATA  ACQUISITION  PROGRAM  FOR 

2  REM  CURRENT  PULSE  LIGHTNING  SIMULATION 

3  REM  AND  SHOCK-EXCITATION  LIGHTNING  SIMULATION 

10  PRINT  'WAVEFORMS  FROM  TAPE  CY  OR  ND';\INPUT  Y$ 

11  IF  Y*-'Y'  THEN  GOT..,  4400 

12  LOAD  'GPI'.'MT' 

13  SIFLIN  BO.'IFC' 

14  SIFCOM  90. 'OCL' 

15  SIFTO  80.500 

20  PRINT  'INITIALIZE  SYSTEM  CY  OR  N)  ';\INPUT  Y$ 

21  IF  Y*-'Y'  THEN  GOTO  1000 

22  GOTO  10S0 

100  REM  START  OF  PROGRAM******************** 

101  GOSUB  4000NREM  SET  UP  FILE  NAME  ANO  NUMBER 

110  PRINT  'SET  UP  DIGITIZERS  CY  OR  N)';\INPUT  Y* 

111  IF  Y«-'Y'  THEN  GOSUB  2000 

115  TT*-FT*\PRINT  'TEST  NAME  ';FT*;\INPUT  FT*\IF  FTS-'S'  THEN  FT*-TT 

$ 

116  PRINT  FT*\IF  FT*-'R'  THEN  GOTO  1050 

117  PRINT  'ARM  SYSTEM  CY  OR  N)':\INPUT  Y* 

118  IF  Y*-'N'  THEN  GOTO  125 

120  FOR  1-33  TO  34 

121  PUT  'ARM  A.B'  INTO  80.1.96 

122  PUT  'POS  0.0'  INTO  a0.I.96*i_ 

123  PUT  'POS  0.0'  INTO  80. I . 96*2~_ 

124  NEXT  I 

125  PRINT  'DIGITIZERS  ARMED  CY  TO  READ)  CN  TO  REARM) ' ; \ INPUT  Y* 

126  IF  Y*-'N'  THEN  GOTO  120 

127  GOSUB  2100 

130  PRINT  'GRAPH  DATA  CY  OR  N)';\INPUT  Y$ 

131  IF  Y*-'Y'  THEN  GOSUB  3000 

135  PRINT  'SAVE  DATA  CY  OR  N)';\INPUT  Y$ 

136  IF  Y*-'Y'  THEN  GOSUB  4100 

137  PRINT  'INDIVIDUAL  PLOTS  CY  OR  N)';\INPUT  YS 

138  IF  YS-'Y'  THEN  GOSUB  3200 

139  GOTO  115 

1000  REM  INITIALIZE  SYSTEM  *********************************** 

1001  SI*-'SE-9' 

1002  SI-VALCSI*) 

1005  LA-32\TA-64\SA-96 

1006  PL$-'RIN  LOW;  CPL  AC;  POS  0.0;  V/D  .1' 

1007  PR*-'RIN  LOW;  CPL  AC;  POS  0.0;  V/D  .1' 

1008  TA$-'TMBS  A;  MOOE  PRE.512»LTC  RIGHT ;  SBPT  O.'iSI* 

1009  TB*-'TMBS  B;MODE  PRE.512;LTC  RIGHT ; SBPT  0. '&.SI* 

1010  TR<-'LTC  RIGHT; CPL  DC; SRC  EXT; LEV  25' 

1020  TT«-OT*\PRINT  'DATE  ';OTi;\INPUT  DT$\IF  DT*-'S'  THEN  DTS-TTS 

1021  TIME  TT*\PRINT  'TIME  ';TT$;\INPUT  TM*\IF  TM*-'S'  THEN  TM*-TT* 
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1022  SETTIM€  TMJ 

1023  SETDATE  OT* 

1024  TIME  TSNOATE  OS 

1025  PRINT  'TODAY  IS  ':DS:'  TIME  IS  ';TS 
1027  APS-'N' 

1030  REM  SET  UP  WAVEFORM  ARRAYS 

1031  DELETE  AA. BB. CC. DD. WA. WB. WC. WD 

1032  HS-'SEC'\VVS-'VOLTS'\VAS-' AMPS' 

1040  WAVEFORM  WA  IS  AAC2047) . SI . HS. VAS 

1041  WAVEFORM  WB  IS  BBC 2047) , SI. HS. VAS 

1042  WAVEFORM  WC  IS  CCC2047) . SI . HS. VAS 

1043  WAVEFORM  WD  IS  DDC2047) . SI . HS. VAS 

1050  REM  NAME  AND  SCALE 

1051  ZAS-NASNPRINT  '#1  NAME  -  'jZAS»  UNPUT  NASMF  NAS-'S'  THEN  NAS-Z 
A* 

1052  PRINT  '#1  SCALE  -'»SAJ\INPUT  ZASUF  ZASo'S'  THEN  SA-VALCZAS) 

1054  ZAS-NBSXPRINT  '#2  NAME  -  ' > ZASj \ INPUT  NBSMF  NBS«'S'  THEN  NBS-Z 
A* 

1055  PRINT  '#2  SCALE  -'sSBi UNPUT  ZASUF  ZASo'S'  THEN  SB-VALCZAS) 

1057  ZAS-NCSXPRINT  '#3  NAME  -  'jZASj UNPUT  NCSUF  NCS-'S'  THEN  NCS-Z 

AS 

1058  PRINT  '#3  SCALE  -'JSCiUNPUT  ZASUF  ZAS<>'S'  THEN  SC-VALCZAS) 

1060  ZAS-NOSXPRINT  '#4  NAME  -  'jZASj UNPUT  NDSXIF  NDS-'S'  THEN  NDS-Z 
AS 

1061  PRINT  '#4  SCALE  -'iSDi UNPUT  ZAS\IF  ZAS<>'S'  THEN  SD-VALCZAS) 

1063  PRINT  CHR (  27  )  J  CHR C  9 1 )  J ' 2 J '  _ 

1064  PRINT  CHR C 27) J  CHRC91 ) J ' 1H' 

1065  PRINT  '#1 ' , NAS, SA 

1066  PRINT  '#2' . NBS. SB 

1067  PRINT  '#3' , NCS. SC 

1068  PRINT  '#4' . NDS. SO 
1099  GOTO  100 

2000  REM  SET-UP  DIGITIZERS  ***************************** 

2002  UT-LA-1 

2004  PUT  PLS  INTO  iO.UT. 96*1 
2006  PUT  PRS  INTO  iO.UT. 96-2 
2008  PUT  TAS  INTO  iO.UT. 96 
2010  PUT  TBS  INTO  80.UT.96 

2012  PUT  TRS  INTO  iO.UT. 96 

2013  SIFCOM  90.UT. 96. 'GTL' 

2016  UT-UT+1 

2018  IF  UT<LA*3  THEN  GOTO  2004 

2019  PRINT  'MAKE  ANY  CHANGES  NOW' 

2020  RETURN 

2100  REM  READ  DATA  AND  SCALE  ********************************* 

2101  TIME  TSXDATE  DS 

2102  UT-1NG0SUB  2150 

2103  AA-A-128 
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2104  VA-VALCSEGCVA».4.LENCVA*))) 

2105  BB-8-12S 

2106  VB-VAL(SEG(VB$. 4. LENC VB$) 3 ) 

2107  UT-2NG0SUB  2150 

2108  CC-A-128 

2109  VC-VALCSEGCVA*.4.LENCVA*))J 

2110  00-B-128 

2111  V0-VALCSEG(VB*.4.LEN(VB*))) 
2118  RETURN 

2150  REM  READ  DIGITIZER  CUT-UNIT 


2155  PRINT  UTXOELETE  A, B 

2160  PUT  'READ  A'  INTO  90.32-KJT.96 

2165  READBI  A  FROM  80. 64+UT. 96 

2166  PUT  ' V/D7'  INTO  90. 32+UT. 96+1 

2167  GET  VAS  FROM  90. 64+UT.  96-^1 
2170  PUT  'READ  B'  INTO  80. 32+UT. 96 

2175  READS  I  B  FROM  80.  64-HJT,  96 

2176  PUT  'V/D7'  INTO  90. 32-HJT,  96-*-2 

2177  GET  VB*  FROM  80. 64-KJT.  96-*-2 
2180  RETURN 

3000  REM  GRAPH  DATA  *****h.«**«*******h.*hh 

3001  PRINT  CHRC27) J  CHRC91 )  J '2J' 

3002  PRINT  CHRC27J i CHRC91 ) s ' 1 s 1H' 

3003  PRINT  CHRC27) J ' 1 ' 

3004  PAGE 

3005  DELETE  GG.WG  _ 

3006  WAVEFORM  WG  IS  GGC2047) . SI . H*. GV« 

3007  UT-1 

3008  GV«-" 

3010  GOTO  UT  OF  3020.3025.3030.3035.3040 

3020  VIEWPORT  100.450.500.700 

3021  SG-CSA-VA J/32 

3022  GG-AA-SG 

3023  GV*-NA* 

3024  GOTO  3100 

3025  VIEWPORT  550.900.500.700 

3026  SG-CSB*VB)/32 

3027  GG-BB-SS 

3028  GVt-NBt 

3029  GOTO  3100 

3030  VIEWPORT  100.450.150.350 

3031  SG-(SC*VC)/32 

3032  GG-CC*SG 

3033  GV8-NC* 

3034  GOTO  3100 

3035  VIEWPORT  550.900.150.350 

3036  SG-CSD*VD)/32 

3037  GG-OD*SG 
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3038  GV»-ND* 

3039  GOTO  3100 

3040  RESETG 

3041  MOVE  10.40NPRINT  'TEST  NAME' 

3042  MOVE  400. 40XPRINT  'FILE  NAME'.'  DATE'.'  TIME' 

3043  MOVE  10.10XPRINT  FT* 

3044  MOVE  400, 10NPRINT  FF*.D*.T* 

3045  IF  AP$-'N'  THEN  GOTO  3050 

3046  PRINT  CHR ( 27 1 J  CHR (231 

3047  GOTO  3051 

3050  WAIT 

3051  PAGE 

3052  PRINT  CHR (271 J '2' 

3053  RETURN 

3100  REM  GRAPH  DATA  ++^+**+*+*^+*+++++*~+~++*****+* 

3102  PRINT  CHR<27) J ' 1 ' 

3105  WINDOW  0, 51*2000. SG*C-12B1 . SG*128 

3106  SETGR  VIEW. WIND, TICS  5.4.GRAT  4, 4. 2, 2 

3110  GRAPH  WG 

3111  UT-UT+1 

3112  GOTO  3010 

3200  REM  GRAPH  INDIVIDUAL  PLOTS 

3201  PRINT  'SELECT  PLOT  (1-2-3-41  CO  TO  RETURN 1 ' » \ I NPUT  GN\IF  GN-0  TH 
EN  RETURN 

3202  PRINT  CHR (271 J  CHR (91 1  * '2J' 

3203  PRINT  CHR ( 27 1 »  CHR (911* ' 1 H '  _ 

3204  PRINT  CHR (271 s ' 1 ' 

3205  PAGE 

3206  DELETE  GG.WG 

3207  WAVEFORM  WG  IS  GGC20471 . SI . H*. GV$ 

3210  GOTO  GN  OF  3215.3220.3225.3230 

3215  SG-(SA*VAl/32 

3216  GG«AA*SG 

3217  GVS-NA* 

3218  GOTO  3234 

3220  SG-(SS*VBl/32 

3221  G G-BB*SG 

3222  GV*-NB* 

3223  GOTO  3234 

3225  SG-(SC*VCl/32 

3226  GG-CC*SG 

3227  GV*-NC* 

3228  GOTO  3234 

3230  SG- (SO* V01/32 

3231  GG-DD*SG 

3232  GV*-NO* 

3234  GG-GG-MEA(GG(0«5011 

3235  VIEWPORT  100.900.150.700 


-112- 


3236  WINDOW  0. SI*2000. SG*(-128) , SG*128 

3237  SETGR  VIEW. WIND, TICS  5.4.  GRAT  4,  4. 2. 2 

3238  GRAPH  WG 

3240  RESETG 

3241  MOVE  10, 40\PRINT  'TEST  NAME' 

3242  MOVE  400.40NPRINT  'FILE  NAME','  DATE','  TIME' 

3243  MOVE  10. 10NPRINT  FTS 

3244  MOVE  400. 10XPRINT  FFS.OS.TS 

3245  WAIT 

3246  PAGE 

3247  PRINT  CHRC27) s '2' 

3250  PRINT  'INT  OR  FFT  OR  N';\INPUT  YS 

3251  IF  Y$-'INT'  THEN  G0SU8  3300 

3252  IF  Y$«'FFT'  THEN  GOSUB  3400 

3253  GOTO  3200 

3300  DELETE  RL.RR 

3301  WAVEFORM  RR  IS  RLC2047) , SI . HRS. VRS 

3302  INT  WG. RR 

3303  HRS**HS 

3304  VRS-GVS&'  INTEGRATED' 

3305  GOTO  3430 

3400  REM  FFT  OF  DISPLAYED  WAVEFORM 
3405  DELETE  RL.IM.RR.II 

3410  WAVEFORM  RR  IS  RLC 1024) . FI . HRS. VR$ 

3415  WAVEFORM  II  IS  IMC 1024) . FI . HIS. VIS 
3420  RFFT  WS.RR.II  — 

3425  POLAR  RR.II 

3426  RLCO)-MIN(RL) 

3427  HRS-'FREQ  IN  MEG  HZ' 

3430  PRINT  CHRC27) *  CHRC91 )  J '2J' 

3435  PRINT  CHRC27) : CHRC91 ) ; ' 1H' 

3440  PRINT  CHRC27) J ' 1 ' 

3445  PAGE 

3450  VIEWPORT  100.300.150,700 
3460  SETGR  VIEW.SRAT  4. 4. 2. 2 
3470  GRAPH  RR 
3475  RESETG 

3480  MOVE  10. 40NPRINT  'TEST  NAME' 

3485  MOVE  400. 40XPRINT  'FILE  NAME'.'  DATE','  TIME' 

3490  MOVE  10. 10NPRINT  FTS 

3495  MOVE  400. 10XPRINT  FFS.DS.TS 

3499  WAIT  1000 

3500  WAIT 
3505  PAGE 

3510  PRINT  CHRC27) J '2' 

3515  RETURN 

4000  REM  SET  UP  FILE  NAME  FOR  SAVE 
4002  REM  FNS-FILE  NAME 
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4004  REM  FN-FILE  NUMBER 

4006  REM  FF4-C0MPLETE  NAME  FOR  FILE  TO  BE  SAVED 

4007  TT*-FN*\PRINT  'FILE  -  ';FF*\INPUT  FN* 

4008  IF  FN*-'S'  THEN  GOTO  4040 

4009  IF  LENCFN*)<7  THEN  GOTO  4011 

4010  PRINT  'FILE  NAME  TOO  LQNG'NGOTO  4007 

4011  SS*-'?' 

4012  FOR  1-1  TO  LENCFN*) 

4014  TC«ASCCSEG<FN*. I. I)) 

4016  IF  TC>47  THEN  IF  TOSS  THEN  GOTO  4022 
4018  IF  TC>64  THEN  IF  TC<91  THEN  GOTO  4022 
4020  GOTO  4024 
4022  SS*-SS*8.CHRCTC:> 

4024  NEXT  I 

4026  FN*«SEGCSS*.2.LEN(SS*):> 

4028  PRINT  'FILE  NUMBER' J \INPUT  FN 

4029  IF  FN<1000  THEN  GOTO  4031 

4030  PRINT  'FILE  NUMBER  TOO  BIG'NGOTO  4028 

4031  SSS-'OO'iSTRCFN) 

4032  FF*-FN*i'. 'iSEGCSS*. LENCSSS1-2. LENCSS*)) 

4034  PRINT  'FILE  NAME  ■' s  FF* 

4035  RETURN 

4040  FN*-TT* 

4041  GOTO  4034 
4100  REM  SAVE  DATA 

4109  PRINT  'FILE  NAME  ■'* FF* 

4110  OPEN  #1  AS  MT«/F.FF$  FOR  WRITE  WITH  4 

4111  WRITE  #1.FF*.FN*.FN.FT*.D*.T* 

4112  WRITE  #1 . NA*. WA. SA, VA 

4113  WRITE  #1 . NB*. WB. SB. VB 

4114  WRITE  #1 . NC*. WC. SC. VC 

4115  WRITE  idl.NDi.WO.SO.VO 

4116  FN-FN-*-l 

4117  SS4-'00'J.STRCFN) 

4118  FFi-FN$i'. 'iSEG(SS$, LEN(SS$)-2. LENCSS$) ) 

4119  CLOSE  #1 

4120  PRINT  'NEXT  FILE  IS  ' s  FF$ 

4121  RETURN 

4400  REM  READ  WAVEFORMS  FROM  MAG  TAPE 
4405  DELETE  AA. BB. CC. DD. WA. WB. WC. WD 
4410  WAVEFORM  WA  IS  AAC2047) . SI . H*. VA* 

4415  WAVEFORM  WB  IS  BBC2047) . SI . H*. VA* 

4420  WAVEFORM  WC  IS  CCC2047) . SI . H*. VA$ 

4425  WAVEFORM  WO  IS  DDC2047) . SI . H*. VA* 

4430  PRINT  'AUTO  PLOT  (Y  OR  N)'i\INPUT  AP* 

4431  RELEASE  'MT' 

4435  LOAD  'MT' 

4440  REWIND  MT « 


4445  GOSUB  4000 

4450  OPEN  #1  AS  MT*/F.FF$  FOR  READ 
4455  READ  #1 . FF$. FN*. FN, FT$. 0$, T$ 

4460  READ  #1 . NA$. WA. SA, VA 
4465  READ  #1 . NB*. *B. SB. VB 
4470  READ  #1 . NO*. WC. SC. VC 
4475  READ  4*1 .  ND$.  WD.  SD.  VD 

4480  CLOSE  #1 

4481  SOSUB  3000 

4482  IF  APS-'Y'  THEN  GOTO  4490 

4485  PRINT  *  INDIVIDUAL  PLOTS  (Y  OR  ND '  s \ INPUT  Y$\IF  Y$-'Y'  THEN  GOSU 
B  3200 

4490  FN“FN-*-l\SS$"'00'  &STRCFN) 

4495  FF$-FN*i*. 'iSEGCSSS . LEN(SS$)-2, LENCSSS) ) 

4496  IF  AP$»' Y '  THEN  GOTO  4450 

4S00  PRINT  'NEXT  FILE  CY  OR  N)  ':FF$\INPUT  Y$ 

4505  IF  Y*-'YJ  THEN  GOTO  4450 
4510  GOTO  444S 
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